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ABSTRACT 

Neutron stars, asteroids, comets, cosmic-dust granules, spacecraft, as well as 
whatever other freely spinning body dissipate energy when they rotate about 
any axis different from principal. We discuss the internal-dissipation-caused re- 
laxation of a freely precessing rotator towards its minimal-energy mode (mode 
that corresponds to the spin about the maximal-inertia axis). We show that this 
simple system contains in itself some quite unexpected physics. While the body 
nutates at some rate, the internal stresses and strains within the body oscillate at 
frequencies both higher and (what is especially surprising) lower than this rate. 
The internal dissipation takes place not so much at the frequency of nutation 
but rather at the second and higher harmonics. In other words, this mechanical 
system provides an example of an extreme non-linerity. Issues like chaos and 
separatrix also come into play. The earlier estimates, that ignored non-linearity, 
considerably underestimated the efficiency of the internal relaxation of wobbling 
asteroids and comets. At the same time, owing to the non-linearlity of inelastic 
relaxation, small-angle nutations can persist for very long time spans. The latter 
circumstance is important for the analysis and interpretation of NEAR's data on 
Eros' rotation state. Regarding the comets, estimates show that the currently 
available angular resolution of spacecraft-based instruments makes it possible to 
observe wobble damping within year- or maybe even month-long spans of time. 
Our review also covers pertinent topics from the cosmic-dust astrophysics; in par- 
ticular, the role played by precession damping in the dust alignment. We show 
that this damping provides coupling of the grain's rotational and vibrational de- 
grees of freedom; this entails occasional flipping of dust grains due to thermal 
fluctuations. During such a flip, grain preserves its angular momentum, but the 
direction of torques arising from H 2 formation reverses. As a result, flipping grain 
will not rotate fast in spite of the action of uncompensated H2 formation torques. 
The grains get "thermally trapped," and their alignment is marginal. Inelastic 
relaxation competes with the nuclear and Barnett relaxations, so we define the 
range of sizes for which the inelastic relaxation dominates. 

Subject headings: 

Euler equations, elliptic functions of Jacobi, celestial mechanics, inelastic relax- 
ation, inelastic dissipation, nonlinear dynamics, interstellar medium (ISM), cos- 
mic dust, comets, asteroids, Eros. 



PACS: 96.50.Gn, 96.30.Ys, 96.35. Cp, 96.35.Fs, 05.45.-a, 45.40.-f, 
95.10.Ce, 95.10.Fh, 45.50.Pk, 98.38.Cp, 98.58.Ca 
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1. INTRODUCTION 
1.1. Motivation 

On the 14-th of February 1967 the Soviet Union launched artificial spacecraft Kosmos 
142, to carry out some ionospheric research. The sputnik had the shape of a cross constituted 
by four 15-meter-long rods. A separate container, shaped as a cylinder with hemispheres on 
its ends, was attached in elastic manner to the cross, in a position orthogonal to its plane. 
This block had dimensions of about 1.6 m x 0.8 m, and carried in itself all the scientific 
equipment. It was connected to the cross frame by a joint, and it turned out that this 
perpendicular position of the container was not secured with a sufficient strength. The mission 
planners wanted the satellite to rotate in the plane of the cross at a rate of 2 revolutions per 
second. At a certain point, when the spacecraft was yet gaining rotation, deformation started. 
The cylindrical container overpowered the locking device in the joint, and bent towards the 
plane of the cross-shaped frame. This phenomenon was addressed by Vasil'ev & Kovtunenko 
(1969) who pointed out that the intensity of the effect depends, among other things, upon 
the angular velocity of rotation of the cross frame. In 22 months after that event, on the 
14-th of December 1968, a similar sputnik, Kosmos 259, was launched. Its rotation rate was 
not so swift: less that one revolution per second. This time no deformations of the spacecraft 
was observed, and the mission succeeded. 

The misadventure of Kosmos 142 resulted from the first principles of mechanics: a freely 
rotating top must end up in the spin state that minimises the kinetic rotational energy, 
for a fixed angular momentum. This spin mode can be achieved by one or both of the 
following means: adjustment of shape or/and alteration of the rotation axis. Since the 
Russian spececraft was easily deformable, it "preferred" the first option. Things would go 
differently if the satellite's construction were more rigid. The latter effect was observed back 
in 1958, when the team operating the first American artificial satellite was surprised by some 
unexpected maneuvres that the spacecraft suddenly began to carry out. The satellite, called 
Explorer I, was a very elongated body with four flexible antennas on it. After launching 
and getting to the orbit, it was set to perform steady rotation about its longest dimension. 
However, the flight operators never managed to keep the spacecraft in the designed spin state: 
Explorer persistently deviated from the simple rotation and went into a wobble, exhibiting 
slowly changing complex spin. Naturally, the rotation state was evolving toward that of 
minimal kinetic energy (the angular momentum being fixed). We would remind that the 
state of rotation about the maximal-inertia axis is the one minimising the kinetic energy, 
while spin about the least-inertia axis corresponds to the maximal energy. Hence, one should 
expect that the body will (through some dissipative processes) get rid of the excessive energy 
and will change the spin axis. 

Another example of unsupported rotator subject to internal dissipation is a cosmic-dust 
granule. Due to various spin-up mechanisms (the main of which is catalytic formation of H2 
moleculae on the granule surface (Purcell 1979)), these particles spend most part of their 
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life in rotation. This circumstance gives birth to a whole sequence of subtle effects, which 
determine alignment of the dust relative to the interstellar magnetic field. This alignment 
can be indirectly observed through measuring the polarisation degree of the starlight passing 
through the dust cloud (Lazarian 2000, 1994). It turns out that theoretical description of 
alignment in based on one's knowledge of the granules' typical rotation state: it is important 
whether the dust particles are, predominantly, in their principal spin states or not (Lazarian 
& Efroimsky 1999; Efroimsky 2002). 

Similar to spacecraft and interstellar grains, a comet or an asteroid in a non-principal 
rotation mode will dissipate energy and will, accordingly, return to the stable spin (Prender- 
gast 1958, Burns & Safronov 1973, Efroimsky & Lazarian 1999). Nevertheless, several objects 
were recently found in excited states of rotation. These are asteroid 4179 Toutatis (Ostro et 
al. 1993, Harris 1994, Ostro et al. 1995, Hudson and Ostro 1995, Scheeres et al. 1998, Ostro 
et al. 1999) and comet P/Halley (Jewitt 1997; Peale & Lissauer 1989; Sagdeev et al. 1989; 
Peale 1991; Wilhelm 1987). Quite possibly, tumbling are also comet 46P/Wirtanen (Samaras- 
inha, Mueller & Belton 1996; Rickman & Jorda 1998), comet 29P/Schwachmann-Wachmann 
1 (Meech et al 1993). The existing observational data on asteroid 1620 Geographos may, too, 
be interpreted in favour of wobble (Prokof'eva et al. 1997; Prokof'eva et al. 1996; Ryabova 
2002). 

The dynamics of a freely rotating body is determined, on the one hand, by the initial 
conditions of the object's formation and by the external factors forcing the body out of its 
principal spin state. On the other hand, it is influenced by the internal dissipation of the 
excessive kinetic energy associated with wobble. Two mechanisms of internal dissipation 
are known. The so-called Barnett dissipation, caused by the periodic remagnetisation, is 
relevant only in the case of cosmic-dust-granule alignment (Lazarian & Draine 1997). The 
other mechanism, called inelastic relaxation, is, too, relevant for mesoscopic grains, and plays 
a primary role in the case of macroscopic bodies. Inelastic relaxation results from alternating 
stresses generated inside a wobbling body by the transversal and centripetal acceleration of 
its parts. The stresses deform the body, and inelastic effects cause energy dissipation. 

The external factors capable of driving a rotator into an excited state are impacts and 
tidal interactions, the latter being of a special relevance for planet-crossers. In the case 
of comets, wobble is largely impelled by jetting. Even gradual outgassing may contribute 
to the effect because a spinning body will start tumbling if it changes its principal axes 
through a partial loss or redistribution thereof. Sometimes the entire asteroid or comet may 
be a wobbling fragment of a progenitor disrupted by a collision (Asphaug & Scheeres 1999, 
Giblin & Farinella 1997, Giblin et al. 1998) or by tidal forces. All these factors, that excite 
rotators, compete with the inelastic dissipation that always tends to return the rotator to 
the minimal-energy state. 

Study of comets' and asteroids' rotation states may provide much information about their 
recent history and internal structure. However, theoretical interpretation of the observational 
data will become possible only after we understand quantitatively how inelastic dissipation 
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affects rotation. The kinetic energy of rotation will decrease at a rate equal to that of energy 
losses in the material. Thus, one should first calculate the elastic energy stored in a tumbling 
body, and then calculate the energy-dissipation rate, using the material quality factor Q. 
This empirical factor is introduced for a phenomenological description of the overall effect of 
the various attenuation mechanisms (Nowick & Berry 1972; Burns 1986, 1977; Knopoff 1963; 
Goldreich & Soter 1965). A comprehensive discussion of the Q-factor of asteroids and of its 
frequency- and temperature-dependence is presented in Efroimsky & Lazarian (2000). 



Our review addresses unsupported rotation of rigid and not-entirely-rigid objects. Stated 
differently, we intend to describe behaviour of unsupported rotators of two sorts: ideal (i.e., 
those that are exempt from internal dissipation) and realistic (i.e., those subject to dissipa- 
tion). While the role of dissipative phenomena in the rotating top has become an issue only 
less than half a century ago, complex rotation of an ideal (absolutely rigid) top has been on 
the scientific agenda since, at least, the mid of XVIII-th century. This problem generated 
some of the major mathematical advances carried out by Jacobi, Poinsot and other eminent 
scholars. However, the founding father of this line of study was Euler whose first notes on 
the topic date back to 1750's. 

Leonhard Euler, the most prolific scientist of all times, will forever retain an aura of 
mistery in the eyes of historians. Very few researchers, if any, shared his power of insight 
and his almost superhuman working ability. His life in science consisted of three major 
periods: the first Russian period (which began in 1730, when young Euler retired from 
the Russian navy for the sake of academic career), the Berlin period (that started in 1741, 
when Euler assumed a high administrative position at the Berlin Academy), and the second 
Russian period (which began in 1765, when major disagreements with King Frederich the 
Second moved Euler to accept an invitation from Empress Catherine the Great, to return to 
St. Petersburg). Each of these three periods in Euler's life was marked by numerous scientific 
achievements in all areas of mathematics known at that epoch. 

One of the fields, that grossly benefitted from Euler's attention during his tenure in 
Berlin, was mechanics of an unsupported top. Euler wrote ca 1760 his pivotal result on the 
topic (Euler 1765), his famous equations: 



where ii i2 ,3 are the eigenvalues of inertia tensor of the body. The tensor is defined through 



p being the position of mass element dm relative to the centre of mass of the body. 
Equations (1-1) are merely a reformulation of the simple fact that the torque equals the rate 



1.2. Complex Rotation of a Rigid Body 



d 
dt 



(J, tti) — (Ij — h) ttj £4 = Tj 



(l-l) 




(1-2) 
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of change of the angular momentum. They express this law in a body frame. Among the 
body frames, there exists one (called principal) wherein the inertia tensor is diagonal. In 
(1-1), 12 ,3 are the angular- velocity components as measured in that, principal, coordinate 
system. Quantities are principal-axes-related components of the total torque acting on 
the body. As ever, we shall assume without loss of generality that J 3 > I 2 > h ■ Hence 
the third axis will always be that of major inertia. 

In the body frame, the period of angular-velocity precession about the principal axis 3 
is: r = 2 Ti /u . Evidently, 

Qi/Qi « r- 1 , Ujh « r^e , (1-3) 

e being a typical value of the relative strain that is several orders less than unity. These 
estimates lead to the inequality U Qi <C Iifli , thereby justifying the commonly used 
approximation to Euler's equations 1 : 

h tli - (Ij - I k ) SI, tt k = n . (1-4) 

Naturally, this elegant system of equations has carried since its birth the name of its author. 
It took scientists some more years to understand that the system deserves its given name for 
one more reason: formulae (1.1) are exactly the Euler-Lagrange equations for the Lagrangian 
of an unsupported rigid body. Proof of this fact demands a certain effort. A brief (but still 
not trivial) derivation offered in 1901 by Poincare can be found in the textbook by Marsden 
(2000). 

The Euler equations simplify considerably when two of three moments of inertia ij are 
equal. This is called dynamic symmetry, to distinguish it from the full geometric symmetry. 
Further on, whenever we refer to symmetric top, we shall imply only the dynamic symmetry, 
not the geometric one. This case was addressed by Euler (1765) himself, and later by Lagrange 
(1813) and Poisson (1813). For prolate symmetric rotators (i.e., when 7 3 = I 2 > h ), in 
the absence of external torques, the solution is simple: 

Qi = const , Q 2 = fij_sinu;t , f2 3 = Q±cosut , (1-5) 

where u = (/1//3 — We see that, from the viewpoint of an observer placed on the 

rotating body, the vector of inertial angular velocity f2 describes a circular cone about the 
minor-inertia axis (1) of the body. So does the angular- momentum vector J . Both Q and 



1 Rigorously speaking, in the case when the approximation (1-3) is not satisfied, not only (1-4) fail but 
even equations (1-1) must be somewhat amended. The problem arises from the ambiguity in the choice of 
the body frame. For example, if we prefer to choose the coordinate system wherein the inertia tensor always 
remains diagonal, then the angular momentum will be different from zero in the body frame (and will be of 
order e). If, though, we choose the coordinates in which the angular momentum vanishes, then the inertia 
matrix will no longer be diagonal. With this choice of the body frame, (1-1) should rather be written down 
not in terms of but in terms of all Iij. We shall not elaborate on this issue in our review. 
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J precess about the least-inertia axis at the same rate to = (h/h — l)^i, though at different 
angles from the axis. 

— * 

In an inertial coordinate system, the angular momentum J will not precess, because 
it must keep constant for a free rotator. Instead, it is the least-inertia axis (1) and the 

— * — * 

angular velocity fl that precess about J , in an inertial observer's opinion. (For brevity, we 
denote each vector by one letter, though we, of course, imply that every vector transforms 
appropriately whenever the coordinate system is changed.) 

In the case of oblate (dynamic) symmetry ( 7 3 > I 2 — I\ ) free precession will be 
expressed, in the body frame, by solution 

Qi = Q± cos out , Q 2 = fi±sinwt , Q 3 = const (1-6) 

— * — * 

where to = (h/h — 1)^3 is the mutual rate of circular precession of fl and J about the 
major inertia axis (3). 

The general case of I3 > I 2 > h is quite involved and demands numerics (see Mitchell 
& Richardson (2001); Richardson & Mitchell (1999), and references therein). Still, in the 
absence of external torques the problem can be solved analytically, and Euler coped with it 
(Euler 1765), though to that end he had to introduce functions similar to what we now call 
elliptic integrals. The solution much simplifies when expressed through the elliptic functions 
of Jacobi sn, cn, dn. These were defined and studied by Karl Jacobi (1829) and used by him 
(Jacobi 1849, 1882) to describe free rotation. Jacobi's functions are generalisations of the 
trigonometric ones, in the following sense: while for symmetric prolate and oblate rotators 
the circular precession is expressed by (1-5) and (1-6) correspondingly, in the general case 
-^3 > I2 > h precession is expressed by very similar formulae that contain Jacobi's finctions 
instead of sin and cos : 

Qi — 7 dn (tut, k 2 ) , f2 2 = (3 sn (cot, k 2 ) , VI3 = a cn (cut, k 2 ) , (1-7) 

for J 2 < 2 I 2 T kin , and 

Qi = 7 cn (uit, k 2 ^ , Q 2 = P sn (&t, k 2 ^ , f2 3 = a dn (ujt, k 2 ^ (1-8) 

for J 2 > 2 I 2 T kin . Here the precession rate cu and the parameters a, (3, (3, 7, 7, Co, k 
and k are certain combinations of ^1,2,3, T kin and J 2 . We see that (1-7) is a generalisation 
of (1-5), while (1-8) is that of (1-6). Solution (1-7) approaches (1-5) in the limit of prolate 
symmetry, (J 3 — — > , while solution (1-8) approaches (1-6) in the limit of oblate 
symmetry, (7 2 — h)/h — > . This situation is illustrated by Figure 1. To understand 
this picture, one should keep in mind that two quantities are conserved for a freely-spinning 
body: the angular momentum 

J 2 = i 2 nj + J 2 fi 2 + J 3 2 fi 2 , (1-9) 

(which is conserved exactly) and the kinetic energy 

T km = l - {hQj + I 2 S% + htti} . (1-10) 
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Fig. 1. — The constant-angular- momentum ellipsoid, in the angular- velocity space. The 
lines on its surface are its intersections with the kinetic-energy ellipsoids corresponding to 
different values of the rotational energy. The quasi-stable pole A is the maximal-energy 
configuration, i.e., the state wherein the body spins about its minimal-inertia axis. The 
stable pole C symbolises the minimal-energy state, i.e., rotation about the maximal-inertia 
axis. The angular-velocity vector describes the constant-energy lines, and at the same time 
slowly shifts from one line to another, approaching pole C. The picture illustrates the case 
of an elongated body: 7 3 ~ 7 2 > A- The trajectories are circular near A and remain (in the 
case of an elongated body) virtually circular almost up to the separatrix. The trajectories 
will regain a circular shape only in the closemost proximity of C. 
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(which is conserved only approximately because of the internal dissipation). Expressions (1-9) 
and (1-10) define ellipsoids in the angular-velocity space (fli, Q2, ^3) • Intersection of these 
gives the trajectory described by the tip of vector f2 in the angular-velocity space. On the 
picture we see the angular-momentum ellipsoid (1-9) with lines marked on its surface. These 
lines are its intersection with ellipsoids (1-10) appropriate to several different values of energy 
Tkin ■ For a fixed value of J 2 , i.e., for a fixed angular-momentum surface (1-9), there exists 
an infinite family of kinetic-energy surfaces (1-10) intersecting with it. The largest surface of 
kinetic energy (corresponding to the maximal value of T kin ) will be an ellipsoid that fully 
encloses our angular-momentum ellipsoid and only touches it in point A and its opposite. 
Similarly, the smallest surface of kinetic energy (corresponding to minimal T^n ) will be an 
ellipsoid escribed by our angular-momentum ellipsoid and only touching it from inside, at 
point C and its opposite. For a fixed J , the maximal and minimal possible values of the 
kinetic energy are achieved during rotations about the minimal-inertia and maximal-inertia 
axes, appropriately. In the case of a non-dissipative torque-free rotation, the tip of vector 

— * 

fl will be describing, on Fig. 1, a curve along which the angular- momentum and energy 
ellipsoids intersect (Lamy & Burns 1972). Solution (1-7) is valid for higher energies, i.e., from 
point A through the separatrix. In astronomy such rotations are called LAM ( = Long- Axis 
Modes). Solution (1-8) works for lower energies, i.e., from the separatrix through point C. 
Such rotations are called SAM ( = Short- Axis Modes). Wherever the trajectories on Fig.l, 
i.e, in the space ( Q\ , Q2 , ^3), are almost circular 2 , the solutions (1-7) and (1-8) may be 
approximated by (1-5) and (1-6), correspondingly. In the limit of an oblate rotator, the 
applicability domain of (1-7) will shrink into a point (or, to be more exact, into two points: 
A and its opposite). Similarly, in the limit of a prolate body, the applicability region of (1-8) 
will shrink into two points: C and its opposite. 



1.3. Realistic Rotators 

The formalism developed by Euler and refined by Jacobi might be a perfect tool for 
description of rotation of asteroids, comets, cosmic-dust granules, spacecraft and whatever 
other unsupported rigid rotators, if not for one circumstance, inner dissipation. Because 
of this circumstance, the Euler- Jacobi theory of precession works only for time spans short 
enough to neglect kinetic-energy losses. 

The necessity of internal dissipation follows from the basic principles of mechanics. A 
freely spinning body of a fixed angular momentum has kinetic energy whose values are con- 
strained to lie within a certain bounded range. Hence, from the physical viewpoint, it is very 
natural for this body to be seeking ways of relaxation. In other words, the body must "do 
its best" to get rid of the excessive kinetic energy, in order to approach the minimal-energy 



2 Be mindful that the trajectory in the space ( Oi , £1 2 , ^3 ) being almost circular does not necessarily 
mean that the precession cone of the major-inertia axis about J is circular or almost circular. 
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configuration. Thence the necessity of some dissipation mechanism. 

Two such mechanisms are known. One is relevant only for mesoscopic rotators, like 
interstellar-dust grains, and therefore plays a certain role in the cosmic-dust alignment. This 
is the Barnett dissipation, a phenomenon called into being by periodic remagnetisation of a 
precessing paramagnetic body (Lazarian & Draine 1997). 

The second mechanism, inelastic dissipation, is, too, relevant for mesoscopic grains 
(Lazarian & Efroimsky 1999), and it plays the decisive role in the macroscopic bodies' relax- 
ation. The effect results from the alternating stresses produced inside a wobbling body by 
the time-dependent acceleration of its parts. The stresses deform the body, and the inelastic 
effects cause dissipation of the rotational energy. 

Dissipation entails relaxation of the precession: the major-inertia axis of the body and 

— * — * 

its angular-velocity vector fl tend to align along the angular momentum J. In other words, 
the precession cone described by O about J will be narrowing until f2 aligns along J 
completely. A simple calculation (Efroimsky 2001, Efroimsky 2000, Efroimsky & Lazarian 
2000, Lazarian & Efroimsky 1999) shows that in this case the major-inertia axis of the body 
will align in the same direction, so that, from the body-frame viewpoint, f2 will eventually 
be pointing along this axis. This configuration will correspond to the minimal kinetic energy, 
the angular momentum being fixed. 

An inertial observer will thus see the unsupported body miraculously changing its rota- 
tion axis. This is exactly what happened in 1958 when, to mission experts' surprise, rotating 
satellite Explorer I changed its rotation axis and went into wobble (Thomson 1961). 

This was, probably, the first example of a practical need for a further development of 
the Eulerian theory of a free top, a development that would address an unsupported top 
with dissipation. However, Chandrasekhar realised this already in mid-50s, after having 
been alerted by Kuiper, and asked a postdoc, Kevin Prendergast, to look into that 3 . The 
most general question was (and still is): how many asteroids in the Solar System can be in 
non-principal (i.e., nutating) spin states, and how can this evidence of the impact frequency 
in the main belt? Prendergast in his paper (1958) implied that it is collisions 4 that drive 
asteroids out of the principal state and make them wobble. An important point made by 
Prendergast (1958) was generation of the second harmonic in a symmetrical oblate rotator: 
if a body is precessing at an angular rate u , then the dissipation is taking place not only at 
this frequency but also at double thereof. Prendergast considered only the deceptively simple 
case of symmetrical rotator, and therefore failed to notice the emergence of harmonics higher 
than the second. Besides, the mathematical treatment of the problem, offered in his paper, 
was erroneous in several other aspects. Nevertheless, the fact that he noticed the second 



3 The authors are grateful to Tom Gehrels for providing these historical facts. 

4 The collisions within the main belt became a popular topic much later, in 90-s. See, for example, 
(dell'Oro, Paolicchi, Cellino, Zappala, Tanga, & Michel 2001). 
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harmonic was by itself an important contribution for which Prendergast should be credited. 
His paper was published much ahead of time and, therefore, was forgotten. Independently 
from Prendergast, Lazarian and Efroimsky (1999) came across the second harmonic some 
40 years later. Generation of the higher harmonics was pointed out only in (Efroimsky 
2000). The reason why the important work by Prendergast was not fully appreciated by his 
contemporaries is that back in 50-s the observational astronomy lacked any reliable data on 
wobbling asteroids. So, Prendergast's paper went almost unnoticed, and his successors had 
to start up from scratch. 

The interest in the asteroidal precession re-emerged in 70-s, after the publication of the 
milestone work by Burns & Safronov (1973) that suggested estimates for the relaxation time, 
based on the decomposition of the deformation pattern into bulge flexing and bending, and 
also on the conjecture that "the centrifugal bulge and its associated strains wobble back and 
forth relative to the body as the rotation axis u moves through the body during a wobble 
period." As turned out later, the latter conjecture was a too strong statement, because the 
inelastic dissipation, for the most part of it, is taking place not near the surface but in the 
depth of the body, i.e., not right under the bulge but deep beneath it. Thus, the bulge 
is much like an iceberg tip. This became clear when the distribution of precession-caused 
stresses was calculated, with improved boundary conditions (Efroimsky & Lazarian 2000), 
(Lazarian & Efroimsky 1999) 5 . Burns & Safronov's treatment neglected the nonlinearity, i.e., 
generation of frequencies higher and lower than the nutation rate. The nonlinearity, in fact, is 
essential. Its neglect leads to a large underestimation of the damping rate, because in many 
spin states a considerable input comes from the harmonics (Efroimsky & Lazarian 2000), 
(Efroimsky 2000). The neglect of nonlinearity leads to up to a two-order underestimate of 
the precession-damping rate. 

In the same year, Peale published an article dealing with inelastic relaxation of nearly 
spherical bodies (Peale 1973), and there he did take the second harmonic into account. 

In 1979 Purcell addressed a similar problem of interstellar-grain precession damping. He 
ignored the harmonics and mishandled the boundary conditions upon stresses: in (Purcell 
1979) the normal stresses had their maximal values on the free surfaces and vanished in the 
centre of the body (instead of being maximal in the centre and vanishing on the surfaces). 
These oversights lead to a several-order underevaluation of the dissipation effectiveness and, 
thereby, of the relaxation rate. 



5 This topic will be discussed in section 4. Our treatment of stresses, demonstrated there, is not math- 
ematically rigorous either. It is a polynomial approximation which satisfies the boundary conditions only 
approximately. From the physical viewpoint, it is not worth further refining that treatment, because the 
slight mishandling of the boundary conditions "spoils" the solution much less than the irregularity and 
inhomogeneity of a the realistic body. 
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1.4. Precession damping 

The dynamics of precession relaxation is described by the angular rate of alignment of 
the maximal- inertia axis (3) along the angular momentum J, i.e., by the decrease in angle 9 
between these. In the case of oblate symmetry (when J 3 > I 2 = h), this angle remains 
adiabatically unchanged over the precession period, which makes d6/dt a perfect measure of 
the damping rate (Efroimsky & Lazarian 2000). However, in the general case of a triaxial 
body angle 9 evolves periodically through the precession cycle. To be more exact, it evolves 
almost periodically, and its value at the end of the cycle is only slightly different from that 
in the beginning of the cycle. The relaxation is taking place through accumulation of these 
slight variations over many periods. This is called adiabatic regime, i.e., regime with two 
different time scales: we have a "fast" process (precession) and a "slow" process (relaxation). 
Under the adiabaticity assumption, one may average 9 , or some function thereof, over the 
precession cycle. Then the damping rate will be described by the evolution of this average. 
Technically, it is convenient to use the average of its squared sine (Efroimsky 2000). One can 
write for a triaxial rotator: 

d < sin 2 9 > d < sin 2 9 > dT kin 
dt dT kin dt 

while for an oblate one the expression will look simpler: 

d9\ ( d9 \ dT^ 



(1-11) 



J (oblate) \dTkin) (oblate) dt 

The derivatives d < sin 2 9 > /dTkin and (d 9 / 'dTkin) ( u a te) appearing in (1-11) and (1-12) 
indicate how the rotational-energy dissipation affects the value of < sin 2 9 > (or simply of 
9 , in the oblate case). These derivatives can be calculated from the equations of motion (see 
Efroimsky & Lazarian (2000) and Efroimsky (2000)). The kinetic-energy decrease, dT k i n /dt , 
is caused by the inelastic dissipation: 

dTk^/dt = < dW/dt > , (1-13) 

W being the energy of the alternating stresses, and < ... > denoting an average over a 
precession cycle. (This averaging is justified within the adiabatic approach. For details see 
section III below.) Finally, in the general case of a triaxial top, the alignment rate will read: 

d < sin 2 9 > d < sin 2 9 > d < W > 
dt dT H n dt 

and for a symmetrical oblate top: 



(1-14) 



d9\ fj^_\ d <W> 

^ J (oblate) V dTkin J ( ob i ate ) dt 

Now we are prepared to set out the strategy of our further work. While calculation of 
d( sin 2 9) /dT^n and (d9 / dTkin) oUate is an easy exercise 6 , our main goal will be to find the 



6 See formula (4-18) below and also formulae (A12 - A13) in Efroimsky 2000. 
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dissipation rate d < W > jdt . This quantity will consist of inputs from the dissipation 
rates at all the frequencies involved in the process, i.e., from the harmonics at which stresses 
oscillate in a body precessing at a given rate u . The stress is a tensorial extension of 
the notion of a pressure or force. Stresses naturally emerge in a spinning body due to the 
centripetal and transversal accelerations of its parts. Due to the precession, these stresses 
contain time- dependent components. If we find a solution to the boundary- value problem 
for alternating stresses, it will enable us to write down explicitly the time- dependent part of 
the elastic energy stored in the wobbling body, and to separate contributions from different 
harmonics: 



W(u n ) being the elastic energy of stresses alternating at frequency uj n . One should know 
each contribution W(u n ), for these will determine the dissipation rate at the appropriate 
frequency, through the frequency- dependent empirical quality factors. The knowledge of 
these factors, along with the averages < W(oo n ) > , will enable us to find the dissipation rates 
at each harmonic. Sum of those will give the entire dissipation rate due to the alternating 
stresses emerging in a precessing body. 



Equation (1-16) implements the most important observation upon which all our study 
rests: generation of harmonics in the stresses inside a precessing rigid body. The harmonics 
emerge because the acceleration of a point inside a precessin body contains centrifugal terms 

— * 

that are quadratic in the angular velocity fl. In the simpliest case of a symmetrical oblate 
body, for example, the body-frame-related components of the angular velocity are given 
in terms of sinut and cosut (see formulae (1-5) - (1-6)). Evidently, squaring of f2 will 
yield terms both with sin out or cos lot and with sin loot or cos 2ut . The stresses produced 
by this acceleration will, too, contain terms with frequency ujt as well as those with the 
harmonic 2ut. In the further sections we shall explain that a triaxial body precessing at 
rate u is subject, in distinction from a symmetrical oblate body, to a superposition of 
stresses oscillating at frequencies u) n = nui, the "base frequency" iQ\ being lower than 
the precession rate uj. The basic idea is that in the general, non-oblate case, the time 
dependence of the acceleration and stresses will be expressed not by trigonometric but by 
elliptic functions whose expansions over the trigonometric functions will generate an infinite 
number of harmonics. In subsection 4.3 we shall explain this in more detail. 

The total dissipation rate will be a sum of the particular rates (Stacey 1992) to be 
calculated empirically. The empirical description of attenuation is based on the quality factor 
Q(u) and on the assumption of attenuation rates at different harmonics being independent 
from one another: 




(1-16) 



n 



1.5. Inelastic dissipation caused by complex rotation 



W = W(u n ) 



n 



E 
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u n W (u n ) 

Q{Un) 
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n 



u n < W(u n ) > 
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Q{uj) being the quality factor of the material, and Wo(cu n ) and < W(u) n ) > being the 
maximal and the average values of the appropriate-to-cj n fraction of elastic energy stored in 
the body. This expression will become more general if we put the quality factor under the 
integral, implying its possible coordinate dependence 7 : 



The above assumption of attenuation rates at different harmonics being mutually independent 
is justified by the extreme smallness of strains (typically, much less than 10~ 6 ) and by the 
frequencies being extremely low (10~ 5 — 1CT 3 Hz). One, thus, may say that the problem is 
highly nonlinear, in that we shall take into account the higher harmonics in the expression 
for stresses. At the same time, the problem remains linear in the sense that we shall neglect 
any nonlinearity stemming from the material properties (in other words, we shall assume 
that the strains are linear function of stresses). We would emphasize, though, that the 
nonlinearity is most essential, i.e., that the harmonics u n come to life unavoidably: no 
matter what the properties of the material are, the harmonics do emerge in the expressions 
for stresses. Moreover, as we shall see, the harmonics interfere with one another due to W 
being quadratic in stresses. Generally, all the infinite amount of multiples of lo\ will emerge. 
The oblate case, where only uj\ and 2uj\ show themselves, is an exception. Another 
exception is the narrow-cone precession of a triaxial rotator studied in Efroimsky (2000): in 
the narrow-cone case, only the first and second modes are relevant (and u\ ss uj). 

Often the overall dissipation rate, and therefore the relaxation rate is determined mostly 
by harmonics rather than by the principal frequency. This fact was discovered only recently 
(Efroimsky & Lazarian 2000, Efroimsky 2000, Lazarian & Efroimsky 1999), and it led to a 
considerable re-evaluation of the effectiveness of the inelastic-dissipation mechanism. In some 
of the preceding publications, its effectiveness had been underestimated by several orders of 
magnitude, and the main reason for this underestimation was neglection of the second and 
higher harmonics. As for the choice of values of the quality factor Q , Prendergast (1958) 
and Burns & Safronov (1973) borrowed the terrestial seismological data for Q. In Efroimsky 
& Lazarian (2000), we argue that these data may be inapplicable to asteroids. 

To calculate the afore mentioned average energies < W(oo n ) > , we use such entities as 
stress and strain. As already mentioned above, the stress is a tensorial generalisation of the 
notion of pressure. The strain tensor is analogous to the stretching of a spring (rendered in 
dimensionless fashion by relating the displacement to the base length). Each tensor compo- 
nent of the stress consists of two inputs, elastic and plastic. The former is related to the strain 
through the elasticity constants of the material; the latter is related to the time-derivative 
of the strain, through the viscosity coefficients. As our analysis is aimed at extremely small 
deformations of cold bodies, the viscosity may well be neglected, and the stress tensor will 
be approximated, to a high accuracy, by its elastic part. Thence, according to Landau & 



'In strongly inhomogcneous nutating bodies attenuation may depend on location. 




(1-18) 
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Lifshitz (1976), the components of the elastic stress tensor <jy are interconnected with those 
of the strain tensor like: 

Tr a ( 1 \ 1 

p, and K being the adiabatic shear and bulk moduli, and Tr standing for the trace of a tensor. 

To simplify the derivation of the stress tensor, the body will be modelled by a rectangular 
prism of dimensions 2a x 2b x 2c where a > b > c. The tensor is symmetrical and is 
defined by 

diGij = paj , (1-20) 

a, being the time- dependent parts of the acceleration components, and p aj being the time- 
dependent parts of the components of the force acting on a unit volume 8 . Besides, the tensor 
<7jj must obey the boundary conditions: its product by normal unit vector, a^rij , must 
vanish on the boundaries of the body (this condition was not fulfilled in Purcell (1979)). 

Solution to the boundary-value problem provides such a distribution of the stresses and 
strains over the body volume that an overwhelming share of dissipation is taking place not 
near the surface but in the depth of the body. For this reason, the prism model gives a 
good approximation to realistic bodies. Still, in further studies it will be good to generalise 
our solution to ellipsoidal shapes. The first step in this direction has been made by Molina, 
Moreno & Martinez-Lopez (2002). 

Equation (1-20) has a simple scalar analogue 9 . Consider a non-rotating homogeneous 
liquid planet of radius R and density p . Let g(r) and P(r) be the free-fall acceleration 
and the self-gravitational pressure at the distance r < R from the centre. (Evidently, 
g(r) = (4/3) nG pr .) Then the analogue to (1-20) will read: 

the expression pg{r) standing for the gravity force acting upon a unit volume, and the 
boundary condition being P(R) = 0. Solving equation (1-21) reveals that the pressure 
has a maximum at the centre of the planet, although the force is greatest at the surface. 
Evidently, the maximal deformations (strains) also will be experienced by the material near 
the centre of the planet. 

In our case, the acceleration a of a point inside the precessing body will be given 
not by the free-fall acceleration g(r) but will be a sum of the centripetal and transversal 
accelerations: fix (f2 x r) + f2 x r , the Coriolis term being negligibly small. Thereby, 



8 Nccdless to say, these acceleration components aj are not to be mixed with a which is the longest 
dimension of the prism. 

9 This example was kindly offered to us by William Newman. 
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the absolute value of a will be proportional to that of r , much like in the above example. In 
distinction from the example, though, the acceleration of a point inside a wobbling top will 
have both a constant and a periodic component, the latter emerging due to the precession. 
For example, in the case of a symmetrical oblate rotator, the precessing components of 
the angular velocity f2 will be proportional to sin ujt and cos out , whence the transversal 
acceleration will contain frequency uj while the centripetal one will contain 2uj. The stresses 
obtained through (1-20) will oscillate at the same frequencies, and so will the strains. As we 
already mentioned, in the case of a non-symmetrical top an infinite amount of harmonics will 
emerge, though these will be obertones not of the precession rate uj but of some different 
"base frequency" uj\ that is less than uj. 

Here follows the expression for the (averaged over a precession period) elastic energy 
stored in a unit volume of the body: 



d < W > 1 1 f / 2 fi 1 



dV 2 < € - > = Tfx \{9K ~ 3 J <{Tra) > + < W > 

{ " l + u -l < ^ > + < <4 > + < °ly > + < °l Z > +2 < <, + < + o\ x > j(l- 

where 2 fx/ (9K) — l/3 = —i//(l+u) ~ —1/5, v being Poisson's ratio (for most solids v pa 1/4). 
Naturally 10 , the total averaged elastic energy is given by the integral over the body's volume: 

<W>= X - J dV ■ cxij €ij , (1-23) 

and it must be expanded into the sum (1-16) of inputs from oscillations of stresses at dif- 
ferent frequencies. Each term (W(uj n )) emerging in that sum will then be plugged into the 
expression (1-17), together with the value of Q appropriate to the overtone uj n . 



22) 



2. NONLINEARITY, CHAOS, SEPARATRIX 

2.1. The Origin of the Nonlinearity 

When (1-18) is inserted into (1-14) and (1-15), one can explicitly see the contributions 
to the entire effect, coming from the principal frequency lo\ and from the harmonics uo n = 
nuj\ . When vector Q describes approximately circular trajectories on Figure 1, the principal 
frequency uj\ virtually coincides with the precession rate uj . This doesn't hold, though, when 
fi get closer to the separatrix: there uj\ becomes lower than the precession rate. The analysis 
of the stress and strain distributions, and the resulting expressions for d < W > /dt written 
down in (Efroimsky & Lazarian 2000) and (Efroimsky 2000) shows that the nonlinearity is 



10 Very naturally indeed, because, for example, a xx e xx dV — (<r xx dy dz)(e xx dx) is a product of the re- 
directed pressure upon the re-directed elongation of the elementary volume dV. 
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essential, in that the generation of harmonics is not a high-order effect but a phenomenon 
playing a key role in the relaxation process. In other words, dissipation associated with the 
harmonics is often of the same order as that at the principal frequency. Near the separatrix 
it may be even higher. 

The nonlinearity emerges due to the simple fact that the acceleration of a point within 
a wobbling object contains centrifugal terms that are quadratic in the angular velocity O . 
In neglect of small terms caused by the body deformation, the acceleration will read: 

a = f2xr + flx(f2xr) . (2-1) 

a being the acceleration in the inertial frame, and r being the position of a point. In 
the simpliest case of oblate symmetry, the body-frame-related components of the angular 
velocity are expressed by (1-5) plugging whereof into (2-1) produces terms containing sinujt 
and cos ujt , as well as those containing sin 2ujt and cos 2ujt . The alternating stresses and 
strains caused by this acceleration are linear functions of a and, thus, will also contain the 
second harmonic, along with the principal frequency. Calculation of the stresses, strains, 
and of the appropriate elastic energy W is then only a matter of some elaborate technique. 
This technique (presented in (Efroimsky & Lazarian 2000) and (Efroimsky 2001)) leads to an 
expression for W , with contributions from uj and 2uj explicitly separated. The nonlinearity 
is essential: in many rotation states the 2uj input in (1-14), (1-15). is of order and even 
exceeds that coming from the principal frequency uj. To explain in brief the reason why 
the nonlinearity is strong, we would mention that while the acceleration and the stresses and 
strains are quadratic in the (precessing) angular velocity fl , the elastic energy is proportional 
to the product of stress and strain tensors. Hence the elastic energy is proportional to the 

— * 

fourth power of fl . 



2.2. The Near-Separatrix Slowing-Down of the Precession 

(Lingering Effect) 

In the general case of a triaxial rotator, precession is described by (1-7) or (1-8). The 
acceleration of a point inside the body (and, therefore, the stresses and strains in the material) 
will, according to (2-1), contain terms quadratic in the Jacobi functions. These functions can 
be decomposed in converging series (the so-called nome expansions) over sin'es and cos'ines of 
nv, n being odd integers for sn(ut, k 2 ) and cn(uit, k 2 ) and even integers for dn(cut, k 2 ) . 
Here v is a frequency lower than the precession rate uj : 

4K(k 2 ) being the mutual period of sn and cn. Near the poles v — > uj, while on approach 
to the separatrix v — > 0. When two such expansions get multiplied by one another, they 
produce a series containing all sorts of products like (sinmz4 sinnz/t) , (cosmz/t cosnz/t), 
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and cross terms. Hence the acceleration, stress and strain contain the entire multitude of 
overtones. Even though the further averaging of W over the precession cycle weeds out much 
of these terms, we are eventually left with all the harmonics on our hands. 

As explained in (Efroimsky 2000), higher-than-second harmonics will bring only high- 
order contributions to the precession-relaxation process when the rotation state is described 
by a point close to poles A or C. Put differently, it is sufficient to take into account only the 
frequencies v « uo and 2v 2uj, insofar as the trajectories on Figure 1 are approximately 
circular (i.e., when (1-6) and (1-7) are well approximated by (1-2) and (1-5)). Near the 
separatrix the situation is drastically different, in that all the harmonics become important. 
We thus transit from the domain of essential nonlinearity into the regime of extreme non- 
linearity, regime where the higher harmonics bring more in the process than v or 1v. We 
are reminded, however, that en route to the separatrix we not just get all the multiples of 
the principal frequency, but we face a change of the principal frequency itself: according to 
(2-2), the principal frequency v will be lower than the precession rate lo\ This regime may 
be called "exotic nonlinearity". 

Without getting bogged down in involved mathematics (to be attended to in subsection 
4.3 below), we would just mention here that in the limit of Q approaching the separatrix the 
dissipation rate will vanish, in the adiabatic approximation. This may be guessed even from 
the fact that in the said limit v — > 0. We thus come to an important conclusion that the 
relaxation rate, being very high at a distance from the separatrix, decreases in its closemost 
vicinity. Can we, though, trust that the relaxation rate completely vanishes on the separatrix? 

— * 

No, because in the limit of fl approaching the separatrix the adiabatic approximation will 
fail. In other words, it will not be legitimate to average the energy dissipation over the 
precession cycle, because near the separatrix the precession rate will not necessarily be faster 
than the relaxation rate. A direct calculation shows that even on the separatrix itself the 
acceleration of a point within the body will remain finite (but will, of course, vanish at the 
unstable middle- inertia pole). The same can be said about stress, strain and the relaxation 
rate. So what we eventually get is not a near-separatrix trap but just lingering: one should 
expect relaxing tops to considerably linger near the separatrix. As for Explorer, it is now 
understandable why it easily went wobbling but did not rush to the minimal-energy spin 
state: it couldn't cross the separatrix so quickly. We would call it "lingering effect". There 
is nothing mysterious in it. The capability of near-intermediate- axis spin states to mimic 
simple rotation was pointed out by Samarasinha, Mueller & Belton (1999) with regard to 
comet Hale-Bopp. A similar setting was considered by Chernous'ko (1968) who studied free 
precession of a tank filled with viscous liquid and proved that in that case the separatrix is 
crossed within a finite time interval 11 . 



11 Such problems have been long known also to mathematicians studying the motion with a non-Hamiltonian 
perturbation: the perturbation wants the system to cross the separatrix, but is not guaranteed to succeed in 
it, because some trajectories converge towards the unstable pole (Neishtadt 1980) 
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3. THE APPLICABILITY REALM OF THE ADIABATIC APPROACH 

In the beginning of subsection 1.2 we explained that Euler's equations (1-1) may be 
written down in their approximate form (1-4) in case the nutation-caused deformations of 
the body are negligibly small. Thus it turns out that in our treatment the same phenomenon 
is neglected in one context and accounted for in another: on the one hand, the very process of 
the inelastic dissipation stems from the precession-inflicted small deformations; on the other 
hand, we neglect these deformations in order to write down (1-4). This approximation (also 
discussed in Lambeck 1988) may be called adiabatic, and it remains acceptable insofar as 
the relaxation is slow against rotation and precession. To cast the adiabatic approximation 
into its exact form, one should first come up with a measure of the relaxation rate. Clearly, 
this should be the time derivative of the angle 9 made by the major-inertia axis 3 and 

— * — * 

the angular momentum J . The axis aligns towards J , so 9 must eventually decrease. Be 
mindful, though, that even in the absence of dissipation, 9 does evolve in time, as can be 
shown from the equations of motion. Fortunately, this evolution is periodic, so one may deal 
with a time derivative of the angle averaged over the precession period. In practice, it turns 
out to be more convenient to deal with the squared sine of 9 (Efroimsky 2000) and to write 
the adiabaticity assertion as: 



uj being the precession rate and < ... > being the average over the precession period. The 
case of an oblate symmetrical top 12 is exceptional, in that 9 remains, when dissipation is 
neglected, constant over a precession cycle. No averaging is needed, and the adiabaticity 
condition simplifies to: 



We would emphasise once again that the distinction between the oblate and triaxial cases, 
distinction resulting in the different forms of the adiabaticity condition, stems from the differ- 
ence in the evolution of 9 in the weak-dissipation limit. The equations of motion of an oblate 
rotator show that, in the said limit, 6 stays virtually unchanged through a precession cycle 
(see section IV below). So the slow decrease of 9 , accumulated over many periods, becomes 
an adequate measure for the relaxation rate. The rate remains slow, compared to the rotation 
and precession, insofar as (3-2) holds. In the general case of a triaxial top the equations of 
motion show that, even in the absence of dissipation, angle 9 periodically evolves, though its 
average over a cycle stays unchanged (virtually unchanged, when dissipation is present but 
weak) 13 . In this case we should measure the relaxation rate by the accumulated, over many 
cycles, change in the average of 9 (or of sin 2 9 ) . Then our assumption about the relaxation 
being slow yields (1-15) 



12 Hereafter oblate symmetry will imply not a geometrical symmetry but only the so-called dynamical 
symmetry: I\ = I 2 - 

13 See formulae (Al) - (A4) in the Appendix to Efroimsky 2000. 




(3-1) 




(3-2) 
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The above conditions (3-1) - (3-2) foreshadow the applicability domain of our further 
analysis. For example, of the two quantities, 

+ + Ilnj = J 2 , (3-3) 

hn\ + i 2 n 2 2 + i 3 n 2 3 = 2T km , (3-4) 

only the former will conserve exactly, while the latter will remain virtually unchanged through 
one cycle and will be gradually changing through many cycles (just like (sin 2 9) ). 



4. SYMMETRIC AND ASYMMETRIC ROTATORS 
4.1. Precession of an Oblate Body. 

An oblate body has moments of inertia that relate as: 

h > I 2 = h = I . (4-1) 

We shall be interested in 9, the rate of the maximum-inertia axis' approach to the direction 
of angular momentum J. To achieve this goal, we shall have to know the rate of energy losses 
caused by the periodic deformation. To calculate this deformation, it will be necessary to 
find the acceleration experienced by a particle located inside the body at a point (x, y, z). 
Note that we address the inertial acceleration, i.e., the one with respect to the inertial frame 
(X,Y,Z), but we express it in terms of coordinates x, y and z of the body frame (1,2,3) 
because eventually we shall have to compute the elastic energy stored in the entire body 
(through integration of the elastic-energy density over the body volume). 

The fast motions (revolution and precession) obey, in the adiabatical approximation, the 
simplified Euler equations (1-14). Their solution, with neglect of the slow relaxation, looks 
(Fowles and Cassiday 1986, Landau and Lifshitz 1976), in the oblate case (4-1): 

fii = Q± cos cut , f2 2 = f2j_ sin cut , f2 3 = const (4-2) 

where 

= f2 sin a , Q 3 = fl cos a , (4-3) 

a being the angle made by the major-inertia axis 3 with f2 . Expressions (4-2) show that 

— * 

in the body frame the angular velocity fl describes a circular cone about the principal axis 
3 at a constant rate 

oo = (h-l)n 3 , h = I 3 /I . (4-4) 

So angle a remains virtually unchanged through a cycle (though in the presence of dis- 
sipation it still may change gradually over many cycles). The precession rate uj is of the 
same order as |f2|, except in the case of h — > 1 or in a very special case of fl and J being 
orthogonal or almost orthogonal to the maximal-inertia axis 3. Hence one may call not only 
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the rotation but also the precession "fast motions" (implying that the relaxation process is a 
slow one). Now, let 9 be the angle between the principal axis 3 and the angular- momentum 
J , so that J 3 = J cos 9 and 

tt 3 = ^ = ^ cos 9 (4-5) 



wherefrom 



h h 



u = (h - 1) ^ cos# . (4-6) 



3 



Since, for an oblate object, 

J = h ei + J 2 ^2 e 2 + J 3 fi 3 e 3 = / (^i e x + fi 2 e 2 ) + I 3 Q 3 e 3 , (4-7) 



the quantity £l± = a/^i + fi 2 is connected with the absolute value of J like: 

Q± = ^ sin 9 = ^- h sin 9 , ft = I 3 /I . (4- 



3 



It ensues from (4-3) that fij_/fi 3 = tan a . On the other hand, (4-5) and (4-8) entail: 
^±/^3 = h tan 9 . Hence, 

tana = h tan# . (4-9) 

We see that angle 9 is almost constant too (though it gradually changes through many 
cycles). We also see from (4-7) that in the body frame the angular- momentum vector J 
describes a circular cone about axis 3 with the same rate uj as Q . An inertial observer, 
though, will insist that it is rather axis 3, as well as the angular velocity Q , that is describing 
circular cones around J. It follows trivially from (4-4) and (4-7) that 

Ifl = J - Iue 3 , (4-10) 

— * — * 

whence it is obvious that, in the inertial frame, both Q and axis 3 are precessing about J at 

• — * — * 

rate J/ 1. (The angular velocity of this precession is e 3 = fl x e 3 = (3/1 — w^) x e 3 = 
(3/1) x e 3 .) Interestingly, the rate uo = (h — 1) Q 3 , at which f2 and J are precessing 
about axis 3 in the body frame, differs considerably from the rate J / 1 at which f2 and axis 
3 are precessing around J in the inertial frame. (In the case of the Earth, J/ 1 m 400 a; 
because h is close to unity.) Remarkably, the inertial-frame-related precession rate is energy- 
independent and, thus, stays unchanged through the relaxation process. This is not the case 
for the body- frame- related rate u which, according to (4-6), gradually changes because so 
does 9. 

As is explained above, we shall be interested in the body-frame-related components fii, 2 ,3 
precessing at rate uj about the principal axis 3. Acceleration of an arbitrary point of the 
body can be expressed in terms of these components through formula 

a = a' + f] x f + 2 fi x v' + x (f] x f) , (4-11) 
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where r, v, a are the position, velocity and acceleration in the inertial frame, and r', v' 
and a' are those in the body frame. Here r = r ' . Mind though that v' and a' do 
not vanish in the body frame. They may be neglected on the same grounds as term 1^ 
in (1-11): precession of a body of dimensions ~ I , with period r , leads to deformation- 
inflicted velocities v' ~ el/r and accelerations a' ps el/r 2 , e being the typical order of 
strains arising in the material. Clearly, for very small e , quantities v' and a' are much 
less than the velocities and accelerations of the body as a whole (that are about l/r and 
l/r 2 , correspondingly). Neglecting these, we get, from (4-11) and (4-2), for the acceleration 
at point (x, y, z) : 

a = ei | - Vt 2 ± x cos 2ut + - Q 2 ^ y sin 2ut + z fl 3 h cos ut | + 

+ e 2 |- Q 2 j_ x sin 2u)t — — Q 2 ± y cos 2ut + z Q± Q 3 h sin out | + 

+ e 3 tt 3 (2 - h) (x cos ujt + y sin cot ) } . (4-12) 

Substitution thereof into (1-20), with the proper boundary conditions imposed, yields, for 
an oblate prism of dimensions 2a x 2a x 2c , a > c, to the following approximate 
expressions: 

pQ 2 pQ 2 
&xx = — r^{x 2 — a 2 ) cos2u;t , a yy = -^(y 2 — a 2 ) cos2u;t , a zz — (4-13) 



°xy = I fi ± (a^ 2 + I/ 2 - 2a 2 ) sin 2wt , (4-14) 



P 



^n ± Vl 3 [h (z 2 - c 2 ) + (2 - h) (x 2 - a 2 ) ] cos cot , (4-15) 

= ^Q L Q 3 [h{z 2 - c 2 ) + (2 - h) {y 2 - a 2 ) } sin ut . (4-16) 

In (4-12) - (4-16) we kept only time-dependent parts, because time-independent parts of 
the acceleration, stresses and strains are irrelevant in the context of dissipation. A detailed 
derivation of (4-12) - (4-16) is presented in (Lazarian & Efroimsky 1999). 

Formulae (4-13) - (4-16) implement the polynomial approximation to the stress tensor. 
This approximation keeps the symmetry and obeys (1-20) with (4-12) plugged into it. The 
boundary condition are satisfied exactly for the diagonal components and only approximately 
for the off-diagonal components. The approximation considerably simplifies calculations and 
entails only minor errors in the numerical factors in (4-22). 

The second overtone emerges, along with the principal frequency u , in the expressions 
for stresses since the centripetal part of the acceleration is quadratic in Q . The kinetic 
energy of an oblate spinning body reads, according to (1-10), (4-3), and (4-9): 



Tkin — 77 [I &± + -^3 ^3 2 ] — 



\ sin 2 9 + -J- cos 2 9 



J 2 , (4-17) 
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wherefrom 

- kin 



, n = —r- {h — 1) sin 9 cos 9 = uo J sin # . (4-18) 
dO I 3 

The latter expression, together with (1-15) and (1-22), leads to: 

<W fdT kin \ l dT kin u . j gin (l] . (j . ? (4 _ 1Q) 



dt \ d9 ) dt 
where 

W = W {uj) + W {2uj) = - u 4&T ~ 2 uj ^Stt ~ (4-20) 

~ " 77 {< ^ H > + 2 < W (M >} , (4-21) 

the quality factor assumed to depend upon the frequency very weakly 14 . In the above formula, 
Wq and Wq u are amplitudes of elastic energies corresponding to the principal mode and the 
second harmonic. Quantities < W 1 ^ > = and < W 2u) > = /2 are the appropriate 

averages. Substitution of (4-13) - (4-16) into (1-22), with further integration over the volume 
and plugging the result into (1-15), will give us the final expression for the alignment rate: 

dO/dt - 3 sin^ 63 (c/a) 4 cot 2 + 20 a 2 gj p 

M/dt ~ " ¥ Sm 6 [1 + (c/a) 2 ] 4 (4 " 22) 

where 

= t (4-23) 

is a typical angular velocity. Deriving (4-22), we took into account that, for an oblate 
2a x 2a x 2c prism (where a > c ), the moment of inertia ^3 and the parameter h read: 

'> = Y"' ,e ■ h = 7 = rrb • (4 " 24) 

Details of derivation of (4-22) are presented in (Lazarian & Efroimsky 1999) 15 . 

Formula (4-22) shows that the major-inertia axis slows down its alignment at small 
residual angles. For 9 — * 0, the derivative becomes proportional to 9, and thus, 9 decreases 
exponentially slowly: 9 = Aexp(— (t), where A and ( are some positive numbers 16 . This 



14 The w-dependence of Q should be taken into account within frequency spans of several orders, but is 
irrelevant for frequencies differing by a factor of two. 

15 Our expression (4-22) presented here differs from the appropriate formula in Lazarian & Efroimsky (1999) 
by a factor of 2, because in Lazarian & Efroimsky (1999) we missed the coefficient 2 connecting Wq with 
W i-). 

16 This resembles the behaviour of a pendulum: if the pendulum is initially given exactly the amount of 
kinetic energy sufficient for the pendulum to move up and to point upwards at the end of its motion, then 
formally it takes an infinite time for the pendulum to stand on end. 
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feature, "exponentially slow finish", (which was also mentioned, with regard to the Chandler 
wobble, in Peale (1973), formula (55)) is natural for a relaxation process, and does not lead 
to an infinite relaxation time if one takes into account the finite resolution of the equipment. 
Below we shall discuss this topic at length. 

Another feature one might expect of (4-22) would be a "slow start": it would be good 
if d9/dt could vanish for 9 — > ir/2 . If this were so, it would mean that at 9 = n/2 (i.e., 
when the major-inertia axis is exactly perpendicular to the angular-momentum vector) the 
body "hesitates" whether to start aligning its maximal-inertia axis along or opposite to the 
angular momentum, and the preferred direction is eventually determined by some stochastic 
influence from the outside, like (say) a collision with a small meteorite. This behaviour is 
the simplest example of the famous spontaneous symmetry breaking, and in this setting it 
is desirable simply for symmetry reasons: 9 = n/2 must be a position of an unstable 
equilibrium 17 . Contrary to these expectations, though, (4-22) leaves dd/dt nonvanishing for 
9 — > 7r/2 , bringing an illusion that the major axis leaves the position 9 = tt/2 at a finite 
rate. This failure of our formula (4-22) comes from the inapplicability of our analysis in the 
closemost vicinity of 9 = 7r/2 . This vicinity simply falls out of the adiabaticity domain of 
our treatment, because ui given by (4-6) vanishes for 9 — > n/2 (then one can no longer 
assume the relaxation to be much slower than the precession rate, and hence, the averaging 
over period becomes illegitimate). This is explained in more detail in the next subsection. 

One more situation, that does not fall under the auspices of our analysis, is when u 
vanishes due to (h — 1) — > 0. This happens when c/a approaches unity. According to (4-22), 
it will appear that dO/dt remains nonvanishing for c/a — > 1, though on physical grounds the 
alignment rate must decay to zero because, for c = a, the body simply lacks a major-inertia 
axis. 

All in all, (4-22) works when 9 is not too close to n/2 and c/a is not too close to unity: 

- 6 < (h - 1) ^ cos9 = ) - ( C { a \l n costf . (4-25) 
h 1 + (c/a) 2 

Knowledge of the alignment rate 9 as a function of the precession-cone half-angle 9 enables 
one not only to write down a typical relaxation time but to calculate the entire dynamics of 
the process. In particular, suppose the observer is capable of measuring the precession- cone 
half-angle 9 with an error 5 . This observer will then compute, by means of (4-22), the 
time needed for the body to change its residual half-angle from 9 to 9 — A9 , for A 9 > 5 . 



17 Imaginc a knife freely rotating about its longest dimension, and let the rotation axis be vertical. This 
rotation mode is unstable, and the knife must eventually come to rotation about its shortest dimension, the 
blade being in the horizontal plane. One cannot say, though, which of the two faces of the blade will look 
upward and which downward. This situation is also illustrated by the pendulum mentioned in the previous 
footnote: when put upside down on its end, the pendulum "hesitates" in what direction to start falling, and 
the choice of direction will be dictated by some infinitesimally weak exterior interaction (like a sound, or 
trembling of the pivot, or an evanescent flow of air). 
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This time will then be compared with the results of his further measurements. Below we 
shall show that such observations will soon become possible for spacecraft. 



4.2. Precession of an Exactly Symmetrical Prolate Body. 

At the first glance, the dynamics of a freely-spinning elongated body obeys the same 
principles as that of an oblate one: the axis of maximal inertia will tend to align itself parallel 
to the angular momentum. If we assume that the body is (dynamically) prolate (i.e., that 
^3 = I2 > h ), it will, once again, be convenient to model it by a prism of dimensions 
2a x 2a x 2c, though this time half-size c will be larger than a . Then all our calculations will 
formally remain in force, up to formula (4-18): since the factor h-1 = [l-(c/a) 2 ]/[l + (c/a) 2 } 
is now negative, the right-hand side in (4-18) will change its sign: 

dT T 2 

(h - 1) sm6 cose = -uJ sine . (4-26) 



de h 

Thereby formula (4-19) will also acquire a "minus" sign in its right-hand side: 

de f dTkin\ dTkin , . • , . 

— = — — — ; — = - (uj J sin e) W (4-27) 
dt \ de J dt v 1 v 1 

Formula (4-21) will remain unaltered. Eventually, by using (4-27) and (4-21), we shall arrive 
to a formula that differs from (4-22) only by a sign, provided e still denotes the angle between 
J and the body-frame axis 3 (parallel to dimension 2c): 

de/dt - 3 sin 3 tf 63 (c/a) 4 cot 2 9 + 20 a 2 flgp 

de/dt - ¥ sme (1 + {c/a)r JQ~ (4-28) 

This looks as if axis 3 tends to stand orthogonal to J, which is natural since axis 3 is now 
not the maximal- inertia but the minimal-inertia axis. 

Alas, all this extrapolation is of marginal practical interest, because even a small differ- 
ence between J 2 and J 3 leads to a considerably different type of rotation. This circum- 
stance was pointed out by Black et al (1999) and was comprehensively discussed by Efroimsky 
(2000). Without getting into excessive mathematical details, here we shall provide a simple 
qualitative explanation of the effect. 



4.3. Triaxial and almost prolate rotators 

Typically, asteroids and comets have elongated shapes, and the above formulae derived 
for oblate bodies make a very crude approximation of the wobble of a realistic triaxial body. 
In the case of a triaxial rotator, with I 3 > h > h , the solution to the Euler equations is 
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expressable in terms of elliptic functions. According to Jacobi (1882) and Legendre (1837), 
it will read, for J 2 < 2 I 2 T kin , as 

fii = 7 dn (cut, k 2 ) , fi 2 = P sn (cut, k 2 ) , Sl 3 = a en (cut, /c 2 ) , (4-29) 

while for J 2 > 2 J 2 it will be: 

Qi = 7 cn [cut, A; 2 ) , Q 2 = (3 sn {wt, /c 2 ) , ^3 = a dn (cut, k 2 ) . (4-30) 

Here the precession rate cu and the parameters a, f3, 7 and /c are some algebraic functions 
°f ^1,2,3, T kin and J 2 . For example, k is expressed by 

k = \ h-h P-2hT k ~ 

" V h - h 2I 3 T km - J 2 

for (4-29), and by 

^ h — h 2I 3 Tkin — J 2 

~ V h - h J 2 - 2/iTfcin 

for (4-30). In the limit of oblate symmetry (when I 2 /h — > 1 ), solution (4-30) approaches 
(4-2), while the applicability region of (4-29) shrinks. Similarly, in the prolate-symmetry 
limit ( (^3 — — > 0) the applicability realm of (4-30) will become infinitesimally small. 
The easiest way of understanding this would be to consider, in the space Q,\ , Q 2 , ^3 , the 
angular- momentum ellipsoid J 2 = If Q 2 + l\ Q\ + 7| Q 2 .A trajectory described by 
the angular-velocity vector f2 in the space Q± , Q, 2 , will be given by a line along which 
this ellipsoid intersects the kinetic-energy ellipsoid 2 T kin = l\ f2 2 + I 2 f2 2 + I 3 fi 2 , as on 
Fig.l. Through the relaxation process, the angular- momentum ellipsoid remains unchanged, 
while the kinetic-energy ellipsoid evolves as the energy dissipates. Thus, the fast process, 
nutation, will be illustrated by the (adiabatically) periodic motion of Q along the line of 
ellipsoids' intersection; the slow process, relaxation, will be illustrated by the gradual shift of 
the moving vector Q from one trajectory to another (Lamy & Burns 1972). On Fig.l, we 
present an angular-momentum ellipsoid for an almost prolate body whose angular momenta 
relate to one another as those of asteroid 433 Eros: 1 x 3 x 3.05 (Black et al. 1999). 

— * 

Suppose the initial energy was so high that Q was moving along some trajectory close to 
the point A on Fig.l. This pole corresponds to rotation of the body about its minor-inertia 

— * — * 

axis. The trajectory described by f2 about A is almost circular and remains so until fl 
approaches the separatrix 18 . This process will be described by solution (4-29). In the vicinity 
of separatrix, trajectories will become noticeably distorted. After the separatrix is crossed, 
librations will begin: Q will be describing not an almost circular cone but an elliptic one. 
This process will be governed by solution (4-30). Eventually, in the closemost vicinity of 



18 As we already mentioned above, this trajectory on Fig.l being almost circular does not necessarily mean 
that the precession cone of the major-inertia axis about J is circular or almost circular. 



(4-31) 



(4-32) 
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pole C, the precession will again become almost circular. (This point, though, will never be 
reached because the alignment of Q along (or opposite to) J has a vanishing rate for small 
residual angles: at the end of the relaxation process the relaxation rate approaches zero, so 
that small-angle nutations can persist for long times.) Parameter k shows how far the tip 
of f2 is from the separatrix on Fig.l: k is zero in poles A and C, and is unity on the 

— * 

separatrix. It is defined by (4-31) when fl is between point A and the separatrix, and by 
(4-32) when f2 is between the separatrix and point C. (For details see (Efroimsky 2000).) 

Now suppose that the rotator is almost prolate, i.e., that point C is very closely embraced 
by the separatrices. If such a rotator gets even slightly disturbed, then, dependent upon the 
particular value of (J 3 — ty/h and upon the intensity of the occational disturbance, the 

— * 

vector fl will be either driven away from point C back to the separatrix, without crossing 
it, or will be forced to "jump" over it. Crossing of the separatrix may be accompanied by 
stochastic flipovers 19 . If we assume that (J 3 — ty/h is infinitesimally small, then the 
separatrix will approach point C infinitesimally close, and the faintest tidal interaction or 
collision will be able to push the vector fl across the separatrix. In other words, an almost 
prolate body, during the most part of its history, will be precessing about its minimal-inertia 
axis. 

If in the early stage of relaxation of an almost prolate (73 pa I2) body the tip of vector 
fl is near point A, then its slow departure away from A is governed by formula (9.22) in 
(Efroimsky 2000): 

d (sin 2 g) 
dt ~~ 

~~ ~~7TT~T ( J 3 " J i) i 1 ~ < sin ^>) {^1 [ 2 < sin ^> " 1 " 

1 7 3 — I 2 h / , . 2 D \ \ 

-5 (1_ (sm e)) 



— UJ 



So 2 -! 1 (I - < sin ^>) + 2ujS * TTTT C 1 " < sin ^>)| • ( 4 - 33 ) 



where 



(2I 3 T km - J 2 ) (J 2 - h) iji HJ _ j ) 7j _ j ) , — — 

(4-34) 

9 is the angle between the angular-momentum vector J and the major-inertia axis 3 ; So,i,2 
are some geometrical factors (So = in the case of 7 2 = ^3), and < ... > symbolises an 



19 The flipovers are unavoidable if dissipation of the kinetic energy through one precession cycle is less than 
a typical energy of an occational interaction (a tidal- force-caused perturbation, for example). 
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average over the precession cycle. For < cos 2 9 > not exceeding m 1/7 , this equation has 
an exponentially decaying solution. For c/a = 0.6 that solution will read: 

At « (-A(0)) x 0.08 ^® . (4-35) 

Comparing this with (4-22), we see that at this stage the relaxation is about 15 times faster 
than in the case of an oblate body. 

During the later stage, when f2 gets close to the separatrix, all the higher harmonics 
will come into play, and our estimate will become invalid. How do the higher harmonics 
emerge? Plugging of (4-29) or (4-30) into (4-12) will give an expression for the acceleration 
of an arbitrary point inside the body. Due to (1-20), that expression will yield formulae for 
the stresses. These formulae will be similar to (4-13 - 4-16), but will contain elliptic functions 
instead of the trigonometric functions. In order to plug these formulae for into (1-22), 
they must first be squared and averaged over the precession cycle. For a rectangular prizm 
2a x 2b x 2c , a direct calculation carried out in (Efroimsky 2000) gives: 

(« 2 XX ) = ^ (I - Q) 2 [3 4 (x 2 - a 2 ) 2 E, , (4-36) 

<ai> = £(S + Q)Wy 2 -& 2 ) 2 S 1 , (4-37) 

= £ (1 - ^W* 2 - c>) 2 S x , (4-38) 

((Tra) 2 ) = ^^{(l-Q)(x 2 -a 2 ) 2 + (S + Q)(y 2 -b 2 ) + (1 - S)(z 2 - c 2 )} 2 E 1 , 

(4-39) 

K) = J {(07 + <*" W - b2 ) + (^7 - au)k 2 )(x 2 - a 2 )} 2 ~ 2 , (4-40) 
(<4> = j{(^ + «7)(^ 2 - c 2 ) + {flu - a^x 2 - a 2 )} 2 E 3 , (4-41) 

K) = J {("0 + ^){z 2 - c 2 ) + (a(3 - co^iy 2 - b 2 )} 2 H 4 , (4-42) 

where Q and S are some combinations of Ii, I2, I3 , defined by formula (2.8) in (Efroimsky 
2000). Factors Si,2,3,4 stand for averaged powers of the elliptic functions: 

Ei = ( (sn 2 (u, k 2 ) - < sn 2 (u, k 2 ) > ) 2 ) = 

= < sn 4 (u, k 2 ) > - < sn 2 (u, k 2 ) > 2 , (4-43) 
E 2 = ( (sn(u, k 2 ) cn(u, k 2 ) — ( sn(u, k 2 ) cn(u, k 2 ))) 2 ) = 
= < sn 2 (u, k 2 ) cn 2 (u, k 2 ) > - < sn(u, k 2 ) cn(u, k 2 ) > 2 , (4-44) 
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S 3 = ( ( cn(u, k 2 ) dn(u, k 2 ) — ( cn(u, k 2 ) dn(u, k 2 ) ) ) 2 ) 



< cn (u, k ) dn (u, k ) > — < cn(w, k ) dn(w, k ) > , 



(4-45) 



S 4 = ( (sn(-u, k 2 ) dn(u, k 2 ) — ( sn(u, k 2 ) dn(u, k 2 )) ) 2 ) = 
= < sn 2 (u, k 2 ) dn 2 (u, k 2 ) > — < sn(u, k 2 ) dn(u, k 2 ) > 2 , 



(4-46) 



where averaging implies: 

<...>=-[ ... du , (4-47) 
T Jo 

t being the mutual period of sn and cn and twice the period of dn : 

r = 4 K(k 2 ) = 4 / (1 - k 2 sm 2 ip)~ 1/2 dip . (4-48) 
Jo 

The origin of expressions (4-43 - 4-46) can be traced from formulae (8.4, 8.6 - 8.13) in 
(Efroimsky 2000). For example, expression (4-11), that gives acceleration of an arbitrary 
point inside the body, contains term sn 2 (ut, k 2 ) . (Indeed, one of the components of the 
angular velocity is proportional to sn(...) , while the centripetal part of the acceleration 
is a quadratic form of the angular- velocity components.) The term sn 2 (ut, k 2 ) in the 
formula for acceleration yields a similar term in the expression for a xx . For this reason 
expression (8.6) in (Efroimsky 2000), that gives the time-dependent part of a xx , contains 
sn 2 (...) — < sn 2 (...) > , wherefrom (4-43) ensues. 

Now imagine that in the formulae (4-36 - 4-42) the elliptic functions are presented by 
their series expansions over sines and cosines (Abramovitz & Stegun 1965): 

sn(ut, k 2 ) = — Y^^(u n t) , (4-49) 

n=l ^ 

27T °° * q n/2 

cn(ut,k 2 ) = — 1 n cos(u n t) , (4-50) 

n=l ^ 

_ oo ** „/ 2 

dn(ut, k ) = — + — J2 cos M , (4-51) 

n=0 ^ 

where 

uj n = nu7r/{2K{k 2 )) , q = exp(-7r K \k' 2 ) / K(k 2 )) , k' 2 = l - k 2 (4-52) 
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and the function K(k 2 ) is the complete elliptic integral of the first kind (see (4-48) or (4- 
57)). A star in the superscript denotes a sum over odd n's only; a double star stands for 
a sum over even n's. Insertion of (4-49-4-51) into (4-36-4-42) will produce, after squaring 
of sn, cn, dn , an infinite amount of terms like sin 2 (uj n t) and cos 2 (uj n t) , along with an 
infinite amount of cross terms. The latter will be removed after averaging over the precession 
period, while the former will survive for all n's and will average to 1/2. Integration over the 
volume will then lead to an expression like (1-16), with an infinite amount of contributions 
(W n ) originating from all cu n 's, n = 1, . . . , oo. This is how an infinite amount of overtones 
comes into play. These overtones are multiples not of precession rate uj but of the "base 
frequency" uj\ = urn / {2K{k 2 )) which is lower than uj . Hence the stresses and strains 
contain not only Fourier components oscillating at frequencies higher than the precession 
rate, but also components oscillating at frequencies lower than ui . This is a very unusual 
and counterintuitive phenomenon. 

The above series (called " nome expansions" ) typically converge very quickly, for q <^ 1 . 
Note, however, that q — > 1 at the separatrix. Indeed, on approach to the separatrix we have: 
k — > 1 , wherefrom K(k 2 ) — > oo ; therefore q — > 1 and ui n — > (see eqn. (4-52)). The 
period of rotation (see (4-48)) becomes infinite. (This is the reason why near-separatrix states 
can mimic the principal one.) 

Paper (Efroimsky 2000), addressed relaxation in the vicinity of poles. This case corre- 
sponds to k <C 1 . For this reason we used, instead of (4-49 - 4-51), trivial approximations 
uj 1 pa uj , sn(ut, k 2 ) pa sin(u;t) , cn(ut, k 2 ) ~ cos(wt) , dn(ut, k 2 ) pal. These approxima- 
tions, along with (4-36 - 4-46) enabled us to assume that the terms a 2 z and a 2 z in (4-6) are 
associated with the principal frequency uj , while < a 2 x > , < a 2 y > , < a 2 z > , < (Tra) 2 > 
and a 2 y are associated with the second harmonic 2 ui . No harmonics higher than second 
appeared in that case. However, if we move away from the poles, parameter k will no longer 
be small (and will be approaching unity as we approach the separatrix). Hence we shall have 
to take into account all terms in (4-49 - 4-51) and, as a result, shall get an infinite amount of 
contributions from all ui n 's in (1-22 - 1-16). Thus we see that the problem is very highly non- 
linear. It is nonlinear even though the properties of the material are assumed linear (strains 
e are linear functions of stresses a ). Retrospectively, the nonlinearity originates because 
the dissipation rate (and, therefore, the relaxation rate) is proportional to the averaged (over 
the cycle) elastic energy stored in the body experiencing precession-caused alternating de- 
formations. The average elastic energy is proportional to < a e > , i.e., to < a 2 > . The 
stresses are proportional to the components of the acceleration, that are quadratic in the 
components of the angular velocity (4-29 - 4-30). All in all, the relaxation rate is a quartic 
form of the angular-velocity components that are expressed by the elliptic functions (4-49 - 
4-51). 

A remarkable fact about this nonlinearity is that it produces oscillations of stresses and 
strains not only at frequencies higher than the precession frequency uj but also at frequencies 
lower than uj. This is evident from formula (4-52): the closer we get to the separatrix (i.e., 
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the closer k 2 gets to unity), the smaller the factor n/(2K) , and the more lower-than-u; 
frequencies emerge. 

A quantitative study of near-separatrix wobble will imply attributing extra factors of 
uj n /Q{uj n ) to each term of the series (1-16) and investigating the behaviour of the resulting 
series (1-17). This study will become the topic of our next paper. Nevertheless, some 
qualitative judgement about the near-separatrix behaviour can be made even at this point. 

For the calculation of the dissipation rate (1-17), the value of the average elastic energy 
< W > given by the sum (1-16) is of no use (unless each of its terms is multiplied by 
uj n /Q{uj n ) and plugged into (1-17)). For this reason, the values of the terms < afj > 
entering (1-22) are of no practical value either; only their expansions obtained by plugging 
(4-49 - 4-51) into (4-36 - 4-46) do matter. Nonetheless, let us evaluate < W > near the 
separatrix. To that end, one has to calculate all < a 2 - > 's by evaluating (4-43 - 4-46). Direct 
integration in (4-43 - 4-47) leads to: 

H: = ^{*._ 1 + ^(2-^-3 (|) 2 } , (4-53) 
= 2 = ^ {2 (* 2 - 1) + I (-2-5P)} , (4-54) 

Ss = 5F {§ f 1 + * 2 > + (* 2 - ')} ' (4 " 55) 

e * = sMf (2 * 2 ~ 1 ) + i 1 - • < 4 - 56 > 

K and E being abbreviations for the complete elliptic integrals of the 1st and 2nd kind: 

K = K{k 2 ) = (1 - k 2 sin 2 -^)" 1/2 dip , 

Jo 

(4-57) 

E = E{k 2 ) = r (1 - k 2 sin 2 ?/;) 1/2 dip . 
Jo 

In the limit of k — > 1 , the expression for K will diverge and all will vanish. Then 
all < <y 2 j > will also become nil, and so will < W > . As all the inputs < W(uo n ) > 
in (1-17) are nonnegative, each of them will vanish too. Hence the relaxation slows down 
near the separatrix. Moreover, it appears to completely halt on it. How trustworthy is this 
conclusion? On the one hand, it might have been guessed simply from looking at (4-48): 
since for k — > 1 the period 4K(k 2 ) diverges (or, stated differently, since the frequencies 
uj n in (4-52) approach zero for each fixed n, then all the averages may vanish). On the other 
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hand, though, the divergence of the period undermines the entire averaging procedure: for 
t — > oo , expression (1-13) becomes pointless. Let us have a look at the expressions for 
the angular-velocity components near the separatrix. According to (Abramovits & Stegun 
1965), these expressions may be expanded into series over small parameter (1 — k 2 ) : 

Qi = 7 dn (cot, k 2 ) = 7 {sech (cot) + 
+ ^ (1 - k 2 ) [sinh(ut) cosh(cot) + cut] sech(cot) tanh(cot)^ + O ((1 - k 2 ) 2 ) , (4-58) 

0,2 = (3 sn (cot, k 2 ) = (3 {tanh (cot) + 
+ i (1 - k 2 ) [smh(cot) cosh(cot) - cot] sech 2 (cot)^ + O ((1 - k 2 ) 2 ) , (4-59) 

fi 3 = a cn (cot, A; 2 ) = a {sech(cot) — (4-60) 

- ^ (1 - k 2 ) [smh(cot) cosh(coi) - cot] sech(cot) tanh(cot)^ + O ((I - k 2 ) 2 ) . (4-61) 

These expansions will remain valid for small k 2 up to the point k 2 = 1 , inclusively. It 
doesn't mean, however, that in these expansions we may take the limit of t — > oo . (This 
difficulty arises because this limit is not necessarily interchangeable with the infinite sum of 
terms in the above expansions.) Fortunately, though, for k 2 = 1 , the limit expressions 

Qi = 7 dn (cot, 1) = 7 sech (cot) , (4-62) 
VL 2 = 13 sn (cot, 1) = tanh(cot) , (4-63) 

f2 3 = a cn (cot, 1) = a sech(cot) (4-64) 

make an exact solution to (1-4). Thence we can see what happens to vector f2 when its 
tip is right on the separatrix. If there were no inelastic dissipation, the tip of vector O 
would be slowing down while moving along the separatrix, and will come to halt at one of the 
middle-inertia homoclinic unstable poles (though it would formally take fl an infinite time 
to get there, because f2i and f2 3 will be approaching zero as ~ exp(— cot) ). When O gets 
sufficiently close to the homoclinic point, the precession will slow down so that an observer 
would get an impression that the body is in a simple-rotation state. In reality, some tiny 
dissipation will still be present even for very slowly evolving f2 . It will be present because 
this slow evolution will cause slow changes in the stresses and strains. The dissipation will 
result in a further decrease of the kinetic energy, that will lead to a change in the value of 
k 2 (which is a function of energy; see (4-31) and (4-32)). A deviation of k 2 away from 
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unity will imply a shift of f2 away from the separatrix towards pole C. So, the separatrix 
eventually will be crossed, and the near-separatrix slowing-down does NOT mean a complete 
halt. 

This phenomenon of near-separatrix slowing-down (that we shall call lingering effect) 
is not new. In a slightly different context, it was mentioned by Chernous'ko (1968) who 
investigated free precession of a tank filled with viscous liquid and proved that, despite the 
apparent trap, the separatrix is crossed within a finite time. Recently, the capability of near- 
intermediate-axis rotational states to mimic simple rotation was pointed out by Samarasinha, 
Mueller & Belton (1999) with regard to comet Hale-Bopp. 

We, thus, see that the near-separatrix dissipational dynamics is very subtle, from the 
mathematical viewpoint. On the one hand, more of the higher overtones of the base frequency 
will become relevant (though the base frequency itself will become lower, approaching zero as 
the angular- velocity vector approaches the separatrix). On the other hand, the separartrix 
will act as a (temporary) trap, and the duration of this lingering is yet to be estimated. 

One should, though, always keep in mind that a relatively weak push can help the 
spinning body to cross the separatrix trap. So, for many rotators (at least, for the smallest 
ones, like cosmic-dust grains) the observational reality near separatrix will be defined not 
so much by the mathematical sophistries but rather by high-order physical effects: the solar 
wind, magnetic field effects, etc... In the case of a macroscopic rotator, a faint tidal interaction 
or a collision with a smaller body may help to cross the separatrix. 



5. APPLICATION TO WOBBLING ASTEROIDS 

5.1. Relaxation Rates of Comets and Asteroids 

Knowledge of the alignment rate 9 as a function of the precession-cone half-angle 9 
enables one not only to write down a typical relaxation time but to calculate the entire 
dynamics of the process. In particular, suppose the observer is capable of measuring the 
precession-cone half-angle 9 with an error 5 . This observer will then compute, by means 
of (4-22), the time needed for the body to change its residual half-angle from 9 to 9 — A9 , 
for A 9 > 5 . This time will then be compared with the results of his further measurements. 
Below we shall show that such measurements will soon become possible for spacecraft-based 
observation means. 

First, let us find a typical relaxation time, i.e., a time span necessary for the major-inertia 
axis to shift considerably toward alignment with J . This time may be defined as: 

f 5 d0 

U = J eo d9jdt ' (5_1) 

9 being the initial half-angle of the precession cone (9q < 7r /2), and 5 being the minimal 
experimentally-recognisable value of 9 . A finite 5 will prevent the "slow-finish" divergency. 
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A particular choice of 0q and 5 will lead to an appropriate numerical factor in the final 
expression for t r . As explained in (Efroimsky 2001), t r is not very sensitive to the choice of 
angle #o , as long as this angle is not too small. This weak dependence upon the initial angle 
is natural since our approach accounts for the divergence at small angles ( "exponentially slow 
finish") and ignores the "slow start". Therefore one can take, for a crude estimate, 

9 = n/2 . (5-2) 

For t r it would give almost the same result as, say, 7r/3 or 7r/4 . A choice of 5 
must be determined exclusively by the accuracy of the observation technique: 5 is such a 
minimally recognizable angle that precession within a cone of half-angle 5 or less cannot 
be detected. Ground-based photometers measure the lightcurve-variation amplitude that is 
approximately proportional to the variation in the cross-sectional area of the wobbling body. 
In such sort of experiments the relative error is around 0.01. In other words, only deviations 
from one revolution to the next exceeding 0. 01 mag may be considered real. This corresponds 
to precession-cone half-angles 5 ~ 10° or larger (Steven Ostro, private communication). 
Ground-based radars have a much sharper resolution and can grasp asteroid-shape details 
as fine as 10 m . This technique may reveal precession at half- angles of about 5 degrees. 
NEAR-type missions potentially may provide an accuracy of 0.01° (Miller et al. 1999). For 
a time being, we would lean to a conservative estimate 

5 = 6° , (5-3) 

though we hope that within the coming years this limit may be reduced by three orders due 
to advances in the spacecraft-borne instruments. 

Remarkably, t r is not particularly sensitive to the half-sizes' ratio c/a either, when 
this ratio is between 0.5 - 0.9 (which is the case for realistic asteroids, comets and many 
spacecraft). Our formulae give: 

Hour result) » (1 " 2) for » (2 - 3) 5 = 12 - 18° ; 



t(our result) « (3 - 4) f OT 9 W 7r/4 J (5-4) 

pa i 'q 



t (our result) « (4 - 5) for ~ 7r/2 

pa J Iq 

(Mind though that, according to (4-25), should not approach n/2 too close.) To compare 
our results with a preceding study, recall that according to Burns & Safronov (1973) 



(5-5) 
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The numerical factor in Burns & Safronov's formula is about 100 for objects of small oblate- 
ness, i.e., for comets and for many asteroids. (For objects of irregular shapes Burns and 
Safronov suggested a factor of about 20 in place of 100.) 

This numerical factor is the only difference between our formula and that of Burns & 
Safronov. This difference, however, is quite considerable: for small residual half-angles 9 , our 
value of the relaxation time is two orders smaller than that predicted by Burns & Safronov. 
For larger residual half-angles, the times differ by a factor of several dozens. We see that 
the effectiveness of the inelastic relaxation was much underestimated by our predecessors. 
There are three reasons for this underestimation. The first reason is that our calculation 
was based on an improved solution to the boundary-value problem for stresses. Expressions 
(4-13) - (4-16) show that an overwhelming share of the deformation (and, therefore, of the 
inelastic dissipation) is taking place in the depth of the body. This is very counterintuitive, 
because on a heuristic level the picture of precession would look like this: a centrifugal 
bulge, with its associated strains, wobbles back and forth relative to the body as Q moves 
through the body during the precession period. This naive illustration would make one 
think that most of the dissipation is taking place in the shallow regions under and around 
the bulge. It turns out that in reality most part of the deformation and dissipation takes 
place deep beneath the bulge (much like in the simple example with the liquid planet, that 
we provided in subsection 1.5. The second, most important, reason for our formulae giving 
smaller values for the relaxation time is that we have taken into account the second harmonic. 
In many rotational states this harmonic turns to be a provider of the main share of the entire 
effect. In the expression (63(c/a) 4 cot 2 9 + 20) that is a part of formula (4-22), the term 
63(c/a) 4 cot 2 9 is due to the principal frequency, while the term 20 is due to the second 
harmonic 20 . For c/a belonging to the realistic interval 0.5 — 0.9 , the second harmonic 
contributes (after integration from 9q through 5 ) a considerable input in the entire effect. 
This input will be of the leading order, provided the initial half-angle 9 is not too small 
(not smaller than about 30°). In the case of a small initial half-angle, the contribution of 
the second mode is irrelevant. Nevertheless it is the small-angle case where the discrepancy 
between our formula and (5-5) becomes maximal. The estimate (5-5) for the characteristic 
time of relaxation was obtained in Burns & Safronov (1973) simply as a reciprocal to their 
estimate for 9 ; it ignores any dependence upon the initial angle, and thus gives too long 
times for small angles. The dependence of the dissipation rate of the values of 9 is the third 
of the reasons for our results being so different from the early estimate (5-5). 

Exploration of this, third, reason may give us an important handle on observation of 
asteroid relaxation. It follows from (4-22) that a small decrease in the precession-cone half- 
angle, — A9 , will be performed during the following period of time: 

A , _ ( M) 24 f 1 + ( c / a ) 2 ] 4 1 fj Q (56) 

" 1 ' 3 63 (c/a) 4 cot 2 9 + 20 sin 3 fl a 2 Q 3 p ' 1 ' 



For calculational details, see Lazarian & Efroimsky (1999). 
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For asteroids composed of solid silicate rock, the density may be assumed p ps 2.5 x 
10 3 kg I vn? , while the product in the numerator should be /iQ ~ 1.5 x 10 13 dyne/ cm 2 = 
1.5 x 10 12 Pa as explained in Efroimsky & Lazarian (2000). Burns & Safronov suggested 
a much higher value of 3 x 10 14 dyne/ cm 2 = 3 x 10 13 Pa , value acceptable within the 
terrestial seismology but, probably, inapplicable to asteroids. 

For asteroids composed of friable materials, Harris (1994) suggests the following values: 
p 2 x 10 3 kg/m 3 and pQ k 5 x 10 12 dyne/cm 2 = 5 x 10 11 Pa . Naturally, this 
value is lower than those appropriate for solid rock (Efroimsky & Lazarian 2000), but in our 
opinion it is still too high for a friable medium. Harris borrowed the aforequoted value from 
preceding studies of Phobos (Yoder 1982). Mind, though, that Phobos may consist not only 
of rubble: it may have a solid component in the centre. In this case, a purely rubble-pile 
asteroid may have a lower pQ than suggested by Harris. Anyway, as a very conservative 
estimate for a rubble-pile asteroid, we shall take the value suggested by Harris. 

As for the geometry, let, for example, 9 = it/ 3 and c/a = 0.6 . Then 

If we measure time At in years, the revolution period T = 2ir/Q in hours, the maximal 
half-size a in kilometers, and 9 in angular degrees (A9 = A9° x 1.75 x 10 -2 ), our formula 
(5-6) will yield: 

At(years) = (_A0°) x 1.31 x 10 -7^3pzl = o.33 , (5-8) 

P a (km) a {km) 

where we accepted Harris' values of pQ = 5 x 10 11 Pa and p = 2 x 10 3 kg/m 3 , and the 
angular resolution of spacecraft-based devices was assumed to be as sharp as \A9\ = 0.01° , 
according to Miller et al. (1999). 



5.2. Parameters of Well-Consolidated Asteroids 

Values of the parameters that appear in the above formulae depend both upon the body 
temperature and upon the wobble frequency. The temperature-, pressure- and frequency- 
caused variations of the density p are tiny and may be neglected. This way, we can 
use the (static) densities appropriate to the room temperature and pressure: p( sdlcate ) ^ 
2500 kg m' 3 and p^bon) w 2 000 kg m' 3 . 

As for the adiabatic shear modulus p , tables of physical quantities would provide its 
values at room temperature and atmospheric pressure, and for quasistatic regimes solely. As 
for the possible frequency-related effects in materials (the so-called ultrasonic attenuation), 
these become noticable only at frequencies higher than 10 8 Hz (see section 17.7 in Nowick 
and Berry 1972). Another fortunate circumstance is that the pressure-dependence of the 
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elastic moduli is known to be weak (Ahrens 1995). Besides, the elastic moduli of solids 
are known to be insensitive to temperature variations, as long as these variations are far 
enough from the melting point. The value of /i may increase by several percent when the 
temperature drops from room temperature to 10 K . Dislocations don't affect the elastic 
moduli either. Solute elements have very little effect on moduli in quantities up to a few 
percent. Besides, the moduli vary linearly with substitutional impurities (in which the atoms 
of the impurity replace those of the hosts). However hydrogen is not like that: it enters the 
interstices between the atoms of the host, and has marginal effect on the modulus. As for 
the role of the possible porosity, the elastic moduli scale as the square of the relative density. 
For porosities up to about 20 % , this is not of much relevance for our estimates. 

According to Ryan and Blevins (1987), for both carbonaceous and silicate rocks one may 
take the shear-modulus value /i w 10 10 Pa . 

Theoretical estimation of the Q— factor for asteroids is difficult. As well known from 
seismology, the Q— factor bears a pronounced dependence upon: the chemical composition, 
graining, frequency, temperature, and confining pressure. It is, above all, a steep function 
of the humidity which presumably affects the interaction between grains. The Q— factor is 
less sensitive to the porosity (unless the latter is very high); but it greatly depends upon 
the amount and structure of cracks, and generally upon the mechanical nature of the aggre- 
gate. Whether comets and asteroids are loose aggregates or solid chunks remains unknown. 
We shall address this question in the next subsection. For now, we shall assume that the 
body is not a loosely connected aggregate but a solid rock (possibly porous but nevertheless 
well-consolidated), and shall try to employ some knowledge available on attenuation in the 
terrestial and lunar crust. 

We are in need of the values of the quality factors for silicate and carbonateous rocks. We 
need these at the temperatures about 150 K, , zero confining pressure, frequencies appropriate 
to asteroid precession (10~ 6 — 10~ 4 Hz), and (presumably!) complete lack of moisture. 

Much data on the behaviour of Q— factors is presented in the seismological literature. 
Almost all of these measurements have been made under high temperatures (from several 
hundred up to 1500 Celsius), high confining pressures (up to dozens of MPa), and unavoidably 
in the presence of humidity. Moreover, the frequencies were typically within the range from 
dozens of kHz up to several MHz. Only a very limited number of measurements have been 
performed at room pressure and temperature, while no experiments at all have been made 
thus far with rocks at low temperatures (dozens of K ). The information about the role of 
humidity is extremely limited. Worst of all, only a few experiments were made with rocks 
at the lowest seismological frequences (10 -3 — 10 _1 Hz), and none at frequencies between 
10~ 6 and 10~ 2 Hz, though some indirectly achieved data are available (see Burns 1977, 
Burns 1986, Lambeck 1980, and references therein). 

It was shown by Tittman et al. (1976) that the Q— factor of about 60 measured under 
ambient conditions on an as-received lunar basalt was progressively increased ultimately to 
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about 3300 as a result of outgassing under hard vacuum. The latter number will be our 
starting point. The measurements were performed by Tittman et al. at 20 kHz frequency, 
room temperature and no confining pressure. How might we estimate the values of Q for 
the lunar basalt, appropriate to the lowest frequencies and temperatures? 

As for the frequency-dependence, it is a long-established fact (e.g., Jackson 1986, Karato 
1998) that Q ~ u a with a around 0.25. This dependence reliably holds for all rocks 
within a remarkably broad band of frequencies: from hundreds of kHz down to 10~ 1 Hz. 
Very limited experimental data are available for frequencies down to 10~ 3 if z, and none below 
this threshold. Keep in mind that a being close to 0.25 holds well only at temperatures of 
several hundred Celsius and higher, while at lower temperatures a typically decreases to 
0.1 and less. 

As regards the temperature-dependence, there is no consensus on this point in the geo- 
logical literature. Some authors (Jackson 1986) use a simple rule: 

Q ~ co a exp(A*/RT) , (5-9) 

A* being the apparent activation energy. A more refined treatment takes into account 
the interconnection between the frequency- and temperature-dependences. Briefly speaking, 
since the quality factor is dimensionless, it must retain this property despite the exponential 
frequency-dependence. This may be achieved only in the case that Q is a function not of 
the frequency per se but of a dimensionless product of the frequency by the typical time of 
defect displacement. The latter exponentially depends upon the activation energy, so that 
the resulting dependence will read (Karato 1998): 

Q ~ [uj exp{A*/RT)f , (5-10) 

where A* may vary from 150 - 200 kJ/mol (for dunite and polycristalline forsterite) up 
to 450 kJ/mol (for olivine). This interconnection between the frequency- and temperature- 
dependences tells us that whenever we lack a pronounced frequency-dependence, the temperature- 
dependence is absent too. It is known, for example (Brennan 1981) that at room temperature 
and pressure, at low frequencies (10~ 3 — 1 Hz) the shear Q— factor is almost frequency- 
independent for granites and (except some specific peak of attenuation, that makes Q increase 
twice) for basalts. It means that within this range of frequencies a is small (like 0.1, or so), 
and Q may be assumed almost temperature-independent too. 

Presumably, the shear Q-factor, reaching several thousand at 20 kHz , descends, in 
accordance with (5-10), to several hundred when the frequency decreases to several Hz. 
Within this band of frequencies, we should use the power a ~ 0.25 , as well known from 
seismology. When we go to lower frequencies (from several Hz to the desirable 10~ 6 — 
10~ 4 Hz), the Q-factor will descend at a slower pace: it will obey (5-10) with a < 0.1 . Low 
values of a at low frequencies are mentioned in Lambeck (1980) and in Lambeck (1988) 21 . 



21 See also Knopoff (1963) where a very slow and smooth frequency-dependence of Q at low frequencies is 
pointed out. 
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The book by Lambeck containes much material on the Q-factor of the Earth. Unfortunately, 
we cannot employ the numbers that he suggests, because in his book the quality factor is 
defined for the Earth as a whole. Physically, there is a considerable difference between the 
Q-factors emerging in different circumstances, like for example, between the effective tidal 
Q-factor 22 and the Q-factor of the Chandler wobble). In regard to the latter, Lambeck 
(1988) refers, on page 552, to Okubo (1982) who suggested that for the Chandler wobble 
50 < Q < 100 . Once again, this is a value for the Earth as a whole, with its viscous layers, 
etc. We cannot afford using these numbers for a fully solid asteroid. 

Brennan (1981) suggests for the shear Q— factor the following values 23 : Q^™^ ~ 

250 , Qf sh^r) ~ ^00 . It would be tempting to borrow these values 24 , if not for one 
circumstance: as is well known, absorption of only several monolayers of a saturant may 
dramatically decrease the quality factor. We have already mentioned this in respect to 
moisture, but the fact is that this holds also for some other saturants 25 . Since the asteroid 
material may be well saturated with hydrogen (and possibly with some other gases), its 
Q-factor may be much affected. 

It may be good to perform experiments, both on carbonaceous and silicaceous rocks, 
at low frequencies and temperatures, and with a variety of combinations of the possible 
saturants. These experiments should give us the values for both shear and bulk quality 
factors. The current lack of experimental data gives us no choice but to start with the value 
3300 obtained by Tittman for thoroughly degassed basalts, and then to use formula (5-10). 
This will give us, at T = 150 K and uj = 10~ 5 Hz : 

- 100 (5-11) 

This value of the shear Q-factor for granites and basalts differs from the one chosen in Burns 
and Safronov (1979) only by a factor of 3. For carbonateous materials Q must be surely 
much less than that of silicates, due to weaker chemical bonds. So for carbonaceous rocks, it 
is for sure that 

Q(carb) < 10Q ? ( 5 _ 12 ) 

though this upper boundary is still too high. 

22 A comprehensive study of the effective tidal Q-factors of the planets was performed by Goldrcich and 
Soter (1966) 

23 Brennan mentions the decrease of Q with humidity but unfortunately does not explain how his specimens 
were dried. 

24 These data were obtained by Brennan for strain amplitudes within the linearity range. Our case is 
exactly of this sort since the typical strain in a tumbling body will be about a j [i w p£l 2 a 2 //i . For the size 
a w 1.5 x 10 4 m and frequency not exceeding 10~ 4 Hz, this strain is less than 1CP 6 which is a critical 
threshold for linearity. 

25 like, for example, ethanol (Clark et al. 1980) 
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Consider asteroid 4179 Toutatis. This is an S-type asteroid analogous to stony irons 
or ordinary chondrites, so the solid-rock value of pQ may be applicable to it: pQ ~ 
1.5 x 10 13 dyne/cm 2 = 1.5 x 10 12 Pa . Its density may be roughly estimated as p = 
2.5 x 10 3 kg/m 3 (Scheeres et al. 1998). Just as in the case of (5-8), let us measure the time 
At in years, the revolution period T in hours (T(hours) — 175), the maximal half-radius a 
in kilometers ( a (km) — 2.2 ), and 9 in angular degrees (|A0°| = 0.01). Then (4-35) will 
yield: 

At {years) « 5.1 x 10~ 2 {h 2 ours) = 5.6 x 10 4 years (5-13) 

a (km) 

Presently, the angular-velocity vector f2 of Toutatis is at the stage of precession about A 
(see Fig.l). However its motion does not obey the restriction (cos 2 9} < 1/7 under which 
(4-35) works well. A laborious calculation based on equations (2.16) and (A4) from Efroimsky 
(2000) and on formulae (1), (2) and (11) from Scheeres et al (1998) shows that in the case of 
Toutatis (cos 2 9) « 2/7. Since the violation is not that bad, one may still use (5-13) as the 
zeroth approximation. Even if it is a two or three order of magnitude overestimate, we still 
see that the chances for experimental observation of Toutatis' relaxation are slim. 

This does not mean, though, that one would not be able to observe asteroid relaxation 
at all. The relaxation rate is sensitive to the parameters of the body (size and density) and 
to its mechanical properties ( pQ ) , but the precession period is certainly the decisive factor. 
Suppose that some asteroid is loosely-connected ( pQ = 5 x 10 12 dyne/ cm 2 = 5 x 10 11 Pa 
and p = 2 x 10 3 kg/m 3 ), has a maximal half-size 17 km, and is precessing with a period of 
30 hours, and is not too close to the separatrix. Then an optical resolution of \A9°\ = 0.01 
degrees will lead to the following time interval during which a 0.01° change of the precession- 
cone half-angle will be measurable: 

T 3 

Af („«,«) « 2.12 x 10~ 2 {h ° urs) = 2 years (5-14) 

a (km) 

which looks most encouraging. In real life, though, it may be hard to observe precession 
relaxation of an asteroid, for one simple reason: too few of them are in the states when 
the relaxation rate is fast enough. Since the relaxation rate is much faster than believed 
previously, most excited rotators have already relaxed towards their principle states and 
are describing very narrow residual cones, too narrow to observe. The rare exceptions are 
asteroids caught in the near-separatrix "trap". These are mimicing the principal state. 

5.3. Loosely-Connected Asteroids 

Above we mentioned that the mechanical structure of the small bodies is still in question. 
At present, most astronomers lean toward the so-called rubble-pile hypothesis, in regard to 
both asteroids and comets. The hypothesis originated in mid-sixties (Opik 1966) and became 
popular in the end of the past century (Burns 1975; Weidenschilling 1981; Asphaug & Benz 
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1994; Harris 1996; Asphaug & Benz 1996; Bottke & Melosh 1996a,b; Richardson, Bottke 
& Love 1998; Bottke 1998, Bottke, Richardson & Love 1998; Bottke, Richardson, Michel & 
Love 1999, Pravec & Harris 2000). 

This hypothesis rests on several arguments the main of which is the following: the large 
fast-rotating asteroids are near the rotational breakup limit for aggregates with no tensile 
strength. This is a strong argument, and one would find difficulty to object to it. Still, we 
would object to two other arguments often used in support this hypothesis. One such dubious 
argument is the low density of asteroid 253 Mathilde (about 1.2 g/cm 3 ) . This low density 
(Veverka et al 1998, Yeomans et al 1998) may be either interpreted in terms of the rubble-pile 
hypothesis (Harris 1998), or be put down to Mathilde being perhaps mineralogically akin to 
low-density carbonaceous chondrites, or be explained by a very high porosity. However, in 
our opinion, the word "porous" is not necessarily a synonim to "rubble-pile", even though in 
the astronomical community they are often used as synonims. In fact, a material may have 
high porosity and, at the same time, be rigid. 

Another popular argument, that we would contest, is the one about crator shapes. 
Many colleagues believe that a rigid body would be shattered into smitherines by collisions; 
therefrom they infer that the asteroids must be soft, i.e., rubble. In our opinion, though, 
a highly porous but still consolidated material may stand very energetic collisions without 
being destroyed, if its porous structure damps the impact. It is know from the construction 
engineering that some materials, initially friable, become relatively rigid after being heated 
up (like, for example, asphalt). They remain porous and may be prone to creep, but they 
are, nevertheless, sufficently rigid and well connected. 

We would also mention that, in our opinion, the sharply-defined craters on the surfaces 
of some asteroids witness against the application of the rubble-pile hypothesis to asteroids. 

For these reasons, we expressed in Efroimsky & Lazarian (2000) our conservative opinion 
on the subject: at least some asteroids are well-connected, though we are uncertain whether 
this is true for all asteroids. This opinion met a cold reaction from the community. However, 
it is supported by the recentmost findings. The monolithic nature of 433 Eros is the most 
important of these (Yeomans 2000). Other include 1998KY26 studied in 1999 by Steven 
Ostro and his team: from the radar and optical observations (Ostro et al 1999), the team 
inferred that this 30-meter-sized body, as well as several other objects, is monolithic 26 . 

Still, despite our conservative attitude toward the rubble-pile hypothesis, we have to 
admit that the main argument used by its proponents (the absence of large fast rotators) 
remains valid, and the question why the large fast-spinning asteroids are near the rotational 
breakup limit for loose aggregates is still awaiting its answer. 



26 We should mention here Vesta as a reliable example of an asteroid being a solid body of a structure 
common for terrestial planets: Hubble images of Vesta have revealed basaltic regions of solidified lava flows, 
as well as a deep impact basin exposing solidified mantle. 



-42 - 



As explained in the preceding subsection, to register relaxation of a solid-rock monolith 
may take thousands of years. However, if the body is loosely connected, the inelastic dissi- 
pation in it will be several orders faster and, appropriately, its relaxation rate will be several 
orders higher. 

5.4. Application to Asteroid 433 Eros in Light of Recent Observations 

As already mentioned above, asteroid 433 Eros is in a spin state that is either principal 
one or very close to it. This differs from the scenario studied in (Black et al 1999). According 
to that scenario, an almost prolate body would be spending most part of its history wobbling 
about the minimal-inertia axis. Such a scenario was suggested because the gap between 
the separatrices embracing pole C on Fig.l is very narrow, for an almost prolate body, and 
therefore, a very weak tidal interaction or impact would push the asteroid's angular velocity 
vector fl across the separatrix, away from pole C. This scenario becomes even more viable 
due to the "lingering effect" described in subsection 2.2, i.e., due to the relative slowing down 
of the relaxation in the closemost vicinity of the separatrix. 

Nevertheless, this scenario has not been followed by Eros. This could have happened 
for one of the following reasons: either the dissipation rate in the asteroid is high enough 
to make Eros well relaxed after the recentmost disruption, or the asteroid simply has not 
experienced impacts or tidal interactions for hundreds of millions of years. 

The latter option is very unlikely: currently Eros is at the stage of leaving the main 
belt; it comes inside the orbit of Mars and approaches that of the Earth. It is then probable 
that Eros during its recent history was disturbed by the tidal forces that drove it out of the 
principal spin state. Hence we are left with the former option, one that complies with our 
theory of precession relaxation. The fact that presently Eros is within less than 0.1 degree 
from its principal spin state means that the precession relaxation is a very fast process, much 
faster than believed previously 27 . 

5.5. Unresolved issues 

Our approach to calculation of the relaxation rate is not without its disadvantages. Some 
of these are of mostly aesthetic nature, but at least one is quite alarming. 

As was emphasised in Section 3, our theory is adiabatic, in that it assumes the presence 
of two different time scales or, stated differently, the superposition of two motions: slow 
and fast. Namely, we assumed that the relaxation rate is much slower than the body-frame- 



27 Note that the complete (or almost complete) relaxation of Eros cannot be put down to the low values 
of the quality factor of a rubble pile, because this time we are dealing with a rigid monolith (Yeomans et al. 
2000). 
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related precession rate u (see formulae (3-1) and (3-2)). This enabled us to conveniently 
substitute the dissipation rate by its average over a precession cycle. The adiabatic assertion 
is not necessarily fulfilled when u itself becomes small. This happens, for example, when 
the dynamical oblateness of an oblate (I 3 > I 2 — I± = I) body is approaching zero: 

(h - 1) -> , h = h/ 1 . (5-15) 

Since in the oblate case u is proportional to the oblateness (see (5.4)), it too will approach 
zero, making our adiabatic calculation inapplicable. This is the reason why one cannot and 
shouldn't compare our results, in the limit of (h — 1) — > , with the results obtained by 
Peale (1973) for an almost-spherical oblate body. 

Another minor issue, that has a lot of mathematics in it but hardly bears any physical 
significance, is our polynomial approximation (4-13 - 4-16 , 4-36 - 4-42) to the stress tensor. 
As explained in Section 4, this approximation keeps the symmetry = <7ji and exactly 
satisfies (1-20) with (4-12) plugged in. The boundary conditions are fulfilled exactly for the 
diagonal components of the tensor and approximately for the off-diagonal elements. In the 
calculation of the relaxation rate, this approximation will result in some numerical factor, 
and it is highly improbable that this factor differs much from unity. 

A more serious difficulty of our theory is that it cannot, without further refinement, 
give a reasonable estimate for the duration of the near-separatrix slowing-down mentioned 
in subsection 2.2. On the one hand, many (formally, infinitely many) overtones of the base 
frequency uj\ come into play near the separatrix; on the other hand, the base frequency 
approaches zero. Thence, it will take some extra work to account for the dissipation associated 
with the stresses oscillating at uj\ and with its lowest overtones. (The dissipation due to the 
stresses at these low frequency cannot be averaged over their periods.) 

There exists, however, one more, primary difficulty of our theory. Even though our 
calculation predicts a much faster relaxation rate than believed previously, it still may fail to 
account for the observed relaxation which seems to be even faster than we expect. According 
to the results obtained by NEAR, the upper limit on non-principal axis rotation is better 
than 0.1 angular degree. How to interpret such a tough observational limit on Eros' residual 
precession-cone width? Our theory does predict very swift relaxation, but it also shows that 
the relaxation slows down near the separatrix and, especially, in the closemost vicinity of 
points A and C on Fig.l. Having arrived to the close vicinity of pole C, the angular- 
velocity vector Q must exponentially slow down its further approach to C . For this reason, 
a body that is monolithic (so that its fiQ is not too low) and whose motion is sometimes 
influenced by tidal or other interactions, must demonstrate to us at least some narrow residual 
precession cone. As already mentioned, for the past million or several millions of years Eros 
has been at the stage of leaving the main belt. It comes inside the Mars orbit and approaches 
the Earth. Most probably, Eros experienced a tidal interaction within the said period of its 
history. Nevertheless it is presently in or extremely close to its principal spin state. The 
absence of a visible residual precession indicates that our theory may still be incomplete. In 
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particular, our Q-factor-based empirical description of attenuation should become the fair 
target for criticisms, because it ignores several important physical effects. 

One such effect is material fatigue. It shows itself whenever a rigid material is subject 
to repetitive load. In the case of a wobbling asteroid or comet, the stresses are tiny, but the 
amount of repetitive cycles, accumulated over years, is huge. At each cycle, the picture of 
emerging stresses is virtually the same. Moreover, beside the periodic stresses, there exists 
a constant component of stress. This may lead to creation of "weak points" in the material, 
points that eventually give birth to cracks or other defects. This may also lead to creep, 
even in very rigid materials. The creep will absorb some of the excessive energy associated 
with precession and will slightly alter the shape of the body. The alteration will be such that 
the spin state becomes closer to the one of minimal energy. It will be achieved through the 
slight change in the direction of the principal axes in the body. If this shape alteration is 
due to the emergence of a considerable crack or displacement, then the subsequent damping 
of precession will be performed by a finite step, not gradually. 

Another potentially relevant phenomenon is the effect that a periodic forcing (such as 
the solar gravity gradient) would have on the evolution and relaxation of the precession 
dynamics. It is possible that this sort of forcing could influence the precessional dynamics of 
the body 28 . 

6. APPLICATION TO WOBBLING COMETS 

6.1. Excitation of the Nucleus' Rotation by Reactive Torques 

According to the widely accepted Whipple's model, comet nuclei are conglomerates of 
ice and dust, also called "dirty snowballs" (Whipple 1951). TV-images of P/Halley and 
P/Borrelly nuclei, obtained by spacecraft Vega- 1,2, GIOTTO and DEEP SPACE- 1, respec- 
tively, allow us to suppose that, generally, cometary nuclei have irregular nonconvex shapes 
with typical sizes R* within the range of 1 to 10 kilometers. Numerically generated images 
of comets P/Halley and P/Borelly are presented on Figure 2 and Figure 3. A typical comet 
rotates slower than most asteroids of a like size. Rotational periods of most comets lie within 
the interval from several hours to several days. 

Solar radiation instigates matter sublimation from the comet-nuclei surfaces; this process 
is especially intensive at heliocentric distances r < 3 AU. Sublimated volatiles (H 2 , CO , 
CO2 , . . .) and emansipated dust form an expanding atmosphere called coma. As a result, 
an Earth-based observer will not only register periodic brightness variations of the comet 
(feature common to spinning comets and asteroids), but will also observe effects specific to 
comets solely. These will be the morphological features of the coma. As an obvious 



We are thankful to Daniel Schcercs who drew our attention to this effect. 
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Fig. 2. — The general view of P/Halley nucleus, based on the numeric-shape model developed 
by Stooke and Abergel (1991) who used the Vega -1,2 and Giotto images. 
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example, we would mention the jets: lengthy narrow flows of matter, initiating on the nucleus' 
surface or close to it. A more subtle phenomenon will be twisting of such jets. As was pointed 
out by Whipple (Whipple 1951), evaporating matter should produce a reactive torque, f , 
that acts on the nucleus and alters its spin state 29 . To calculate it, one can use formula 

N 

? = -^<^(R,-xv,-), (6-1) 

where N is the number of faces of the polyhedron approximating the nucleus' shape, Qj is 
the rate of mass ejection by the j-th face, R-,- is the radius vector of the face's centre in the 
body's principal frame of reference, and v,- is the effective velocity of the ejected matter. 
The process of mass ejection depends on local illumination conditions and on the current 
heliocentric distance. Its accurate description is a subject of a rather involved research (Crifo 
& Rodionov 1999). Various empirical models have been employed to study rotation evolution 
of cometary nuclei. For example, to model the relative dependence of the mass ejection rate 
on the helicentric distance, the following expression was suggested in (Marsden, Sekanina & 
Yeomans 1973) and has been widely used since then: 
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Here r = 2.808 AU, Qq is a normalising multiplier (whose its value can be chosen in such a 
way that g(q) = 1, where q stands for the perihelion distance). 

In most publications on the topic, spin evolution of a comet nuclei due to reactive torques 
was studied through numeric integration of the equations of rotation (Wilhelm 1987, Julian 
1990, Peale & Lissauer 1989, Samarasinha & Belton 1995, Jorda & Licandro 2002). A more 
systematic approach to the problem was devoloped in (Neishtadt, Scheeres, Sidorenko & 
Vasiliev 2002), where possible secular phenomena were studied by means of an accurate aver- 
aging procedure (Arnold 1978), and relevant physical parameters, that control the evolution 
of cometary rotation, were extracted. 

Numeric and analytic investigations have revealed certain patterns in the spin evolution. 
For nuclei with fixed active regions, the direction of the angular momentum vector L typically 
spirals either toward the orbital direction of peak outgassing or opposite to it (Neishtadt, 
Scheeres, Sidorenko & Vasiliev 2002; Samarasinha & Belton 1995). The timescale for the 
reactive torque to change the nucleus spin state can be estimated as 



29 There are analogies between the spin-up of comets due to matter ejection and spin-up of grains due to 
ejection of H 2 atoms at the catalytic sites on grain surface. The theory of grain rotation due to such torques 
is discussed in Purccll (1979), Spitzer & McGlynn (1979), Lazarian & Draine (1997). 
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Fig. 3. — P/Borelly nucleus: rough reconstruction based on Deep Space - 1 images. According 
to (Neishtadt, Scheeres, Sidorenko, Stooke & Vasiliev 2002), this nucleus can be approximated 
by a combination of two ellipsoids with major semiaxes 1.6 km, 1.8 km, 3.0 km, and 0.96 km, 
1.08 km, 1.8 km, the distance between their centres being 3.7 km. 
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being the nucleus angular momentum, Q* being the total mass ejection rate at perihelion, 
and V* being the effective exhaust velocity. The product (D*R*) can be interpreted as 
the torque's arm. Here the dimensionless parameter characterises the dependence of 
reactive torque on active zones' distribution over the nucleus surface, while the dimensionless 
parameter is a certain integral characteristic of the nucleus' heating on its orbit: 



Here T and e are comet's orbital period and eccentricity, while v is the nucleus true anomaly. 

Although the long-term history of the nucleus spin state is dominated by reactive torques, 
we should not forget that nuclei are not perfectly rigid bodies. Splitting of comets' nuclei by 
tidal stress, phenomenon that has been observed many a time by astronomers, has inspired 
some theoreticians to speculate that nuclei consist of weakly bonded components. This so- 
called "rubble-pile hypothesis" has been extended by some authors even to asteroids. As 
we already mentioned in one of the preceding sections of this review, this hypothesis has 
its pro's and con's and, therefore, both proponents and critics. What is certain thus far, is 
that at least some comets are indeed weak. For example, two unequal components of comet 
P/Borrelly are connected with an intermediate zone covered by a complicated network of 
cracks. This structure, well visible on the available images, has lead to a hypothesis that the 
P/Borrelly nucleus is assembled of two pieces in loose contact. 

When the comet is not too far from the Sun, time tn is, typically, much shorter than 
the internal-relaxation time. Nevertheless, it is possible that some qualitative properties of 
the nucleus' rotation are associated with non-rigidity. In particular, since the timescale for 
wobble excitation by r is often shorter than the damping time, it is reasonable to suppose 
that many cometary nuclei should be in excited rotational states after passing their perihelia 
(and sometimes even before approaching the perihelia). So far such an excited state has been 
reliably established for P/Halley only 30 . However, as we mentioned above there are several 
more comets whose spin state may be complex. 

The small amount of tumbling comets observed thus far may be put down to one or 
both of the following reasons. The first is the still insufficient precision of our observation 
techniques. The second reason is the possibility of relaxation rate being much higher than 
believed previously. 



30 According to (Bclton, Julian, Anderson & Mueller 1991), the long axis of P/Halley nucleus is inclined 
to the angular-momentum vector L by « 66° and is rotating about L with a period of 3.69 days. The spin 
component about the long axis has a period of 7.1 days. 




(6-4) 
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6.2. Parameters of the Comets 

For the dirty ice Peale and Lissauer (1989) suggest, for Halley's comet, ji ~ 10 10 dyne /cm 2 
= 10 9 Pa while Q < 100 . The latter number, 100, characterises attenuation in the solid 
ice. The question arises whether it is applicable to comets. In the preceding section we men- 
tioned, in regard to asteroids, the "rubble-pile" hypothesis. Despite that hypothesis being 
fortified by various theoretical explanations, our attitude to it still remains conservative. 

In the case of comets, though, this hypothesis may work. It is a well-established fact that 
comets sometimes get shattered by tidal forces. For example, Shoemaker-Levy 9 broke into 21 
pieces on the perijove preceding the impact (Marsden 1993). In 1886 Comet Brooks 2 was rent 
into pieces by the jovian gravity (Sekanina 1982). Comet West disintegrated in 1976 (Melosh 
and Schenk 1993). It seems that at least some comets are loosely consolidated aggregates, 
though we are unsure if this is the case for all comets. Hopefully, our understanding of the 
subject will improve after the Deep Impact mission reaches its goal. 

Loosely-connected comets, presumably, consist not so much of solid ice but more of firn 
(coarse-grained snow), material whose quality factor is of order unity. Hence the estimate 

Qicomet) ^ j (g_ 5 ) 

is, in our opinion, more realistic. So, let us assume fiQ ~ 10 10 dyne /cm 2 = 10 9 Pa. As for the 
average density of a comet, most probably it does not deviate much from 1.5 x 10 3 kg/m 3 . 
Indeed, on the one hand, the major part of the material may have density close to that 
of firn, but on the other hand a typical comet will carry a lot of crust and dust on and 
inside itself. Now, consider a comet of a maximal half-size 7.5 km (like that of Halley comet 
(Houpis and Gombosi 1986)) precessing with a period of 3.7 days « 89 hours (just as Halley 
does 31 ). If we once again assume the angular resolution of the spacecraft-based equipment 
to be |A0°| = 0.01, it will lead us to the following damping time: 

T 3 

&t iyears) « 5.65 x 10~ 5 {h ° urs) = 0.7 year . (6-6) 

a {km) 

This means that the cometary-relaxation damping may be measurable. It also follows from 
(6-6) that, to maintain the observed tumbling state of the Comet P/Halley, its jet activity 
should be sufficiently high. 

6.3. Missions to Comets, and the Possible Nuclei-Dynamics Studies 

One of the main results of our research is that, according to (6-6), wobble damping 
of comets may be registered within a one-year time span or so, provided the best available 



Bclton ct al 1991, Samarasinha and A'Hearn 1991, Peale 1992 
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spacecraft-based devices are used and the comet's spin state is not too close to the separatrix. 
In the coming years a new generation of comet-aimed missions will increase the amount 
and quality of observations. Forthcoming enlargement of the comet database will result in 
deepening of our understanding of comet nuclei rotational dynamics. 

The Stardust mission, that is to visit comet P/Wild 2, will be a fly-by one, and thus will 
not be able to trace the spin-state evolution. Spececraft Deep Impact will approach comet 
9P Tempel 1 on the 3 of July 2005, and will shoot a 500-kg impactor at 10 km/s speed, to 
blast a crater into the nucleus, to reveal its interior. Sadly, though, this encounter too will 
be short. 

Rosetta orbiter is designed to approach 46/P Wirtanen and to escort it for about 2.5 
years (Hubert & Schwehm 1991). Wirtanen seems to be a wobbling comet (Samarasinha, 
Mueller & Belton 1996), (Rickman & Jorda 1998), and one of the planned experiments is 
observation of its rotation state, to be carried out by the OSIRIS camera (Thomas et al 
1998). The mission will start about 1.5 year before the perihelion, but the spin state will be 
observed only once, about a year before the perihelion, at a heliocentric distance of 3.5 AU. 
Three months later the comet should come within 3 AU which is to be a crucial threshold 
for its jetting activity: at this distance outgassing of water will begin. The strongest non- 
gravitational torques emerge while the comet is within this distance from the Sun. It is 
predominantly during this period that wobble is instigated. Hence, it may be good to expand 
the schedule of the spin-state observations: along with the measurement currently planned for 
3.5 AU before perihelion, another measurement, at a similar heliocentrical distance after the 
perihelion, would be useful. The comparison of these observations will provide information 
of the precession increase during the time spent by the comet within the close proximity to 
the Sun. Unfortunately, the Rosetta programme will be over soon afterwards, and a third 
observation at a larger distance will be impossible. The third observation performed well 
outside the 3 AU region might reveal the wobble-damping rate, and thereby provide valuable 
information about the composition and inner structure of the comet nucleus. It would be 
most desirable to perform such a three-step observation by the future escort missions to 
comets. What are the chances of success of such an experiment? On the one hand, the 
torques will not be fully eliminated after the comet leaves the 3 AU proximity of the Sun; 
though the outgassing of water will cease, some faint sublimation of more volatile species 
(like CO, C0 2 , CH 3 OH) will persist for long. On the other hand, the damping rate is high 
enough and a comparison of two observations separated by about a year will have a very 
good chance of registering wobble relaxation. 



7. APPLICATION TO COSMIC-DUST ALIGNMENT 

Inelastic relaxation is also important for cosmic-dust grains of sizes within the interval 
from 3 x 10~ 7 cm through 10~ 4 cm. Such particles get their long axes aligned with respect 
to the interstellar magnetic field (see review by Lazarian 2000 for description of alignment 
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mechanisms). Absorption and emission of light by such particles is charachterised by appro- 
priate effective cross sections which depend upon the direction from which the observer looks 
at the particle. The difference between cross sections appropriate to the directions parallel 
to the long and short axes of the granule yields polarisation of the observed starlight, in case 
this light passes through clouds of aligned dust. Dust particles tend to align their orienta- 
tion relative to the interstellar magnetic field, and it turns out that the polarisation of light 
depends upon the orientation and magnitude of the field. This circumstance gives birth to a 
very powerful technique for magnetic field studies, technique that rests on our understanding 
of the grain alignment, and on our ability to observe this alignment through measuring the 
polarisation of starlight passing through dust clouds. 

Dust interact with the ambient media and rotate fast. There are a number of processes 
that allow grains to spin at a rate much larger than the typical rate of thermal rotation 
Wth ~ (kT gas / Jgrain) 1 ! 2 1 where J grain is a typical angular momentum of the grain. For instance, 
formation of H2 molecules over catalytic sites on grain surface gives rise to the so-called Purcell 
rockets (Purcell 1979), which can spin up grains by ejecting nascent H2 molecules. Under 
such uncompensated torques the grain will spin-up to velocities much higher than u t h, and 
Purcell termed those velocities "suprathermal." 

It is obvious that grain wobble limits the degree to which the grain can be aligned. The 
ideas on the role of this wobble have been evolving through years. Initially it was assumed 
that the distribution of grain rotation among the axes is thermal Iiuf ~ kT gra i n (see Jones 
& Spitzer 1967). Later Purcell (1979) realised that internal relaxation tends to bring grain 

— * 

angular momentum J parallel to the axis of maximal inertia. He discussed both inelastic 
relaxation and a more subtle mechanism that he termed Barnett relaxation. 

The Barnett effect is converse of the Einstein-Haas effect. The Einstein-Haas effect 
is rotation acquired by a paramagnetic body subject to remagnetisation. This happens 
because the flipping electrons transfer the angular momentum (associated with their spins) 
to the lattice. In the case of Barnett effect, the rotating body shares a part of its angular 
momentum with the electron subsystem; this sharing entails magnetisation. A typical value of 
the Barnett-induced magnetic moment is M ~ 10~ 19 c<j(5) erg gauss -1 (where uj^ = cj/10 5 s -1 ). 
Purcell offered the following illustrative explanation of the phenomenon. If a rotating body 
contains equal amount of spin-up and spin-down unpaired electrons, its magnetisation is nil. 
Its kinetic energy would decrease, with the total angular momentum remaining unaltered, if 
some share of the entire angular momentum could be transferred to the spins by turning some 
of the unpaired spins over (and, thus, by dissipating some energy). This potential possibility 
is brought to pass through the said coupling. Another way to understand the effect would be 
to recall that in the body frame of a free rotator the angular-momentum and angular- velocity 
vectors precess relative to the principal axes. Therefore, the Barnett magnetisation will, too, 
precess in the body frame. The resulting remagnetisation will be accompanied with energy 
dissipation (called paramagnetic relaxation). 

An immediate outcome from granule magnetisation is the subsequent coupling of the 
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magnetic moment M (directed along the grain angular velocity) with the interstellar mag- 
netic field B : the magnetic moment precesses about the magnetic line. What is important, 
is that this precession goes at an intermediate rate. On the one hand, it is slower than 
the grain's spin about its instantaneous rotation axis. On the other hand, the precession 
period is much shorter than the typical time scale at which the relative-to-B alignment gets 
established 32 . The latter was proven by Dolginov & Mytrofanov (1976), for magnetisation 
resulting from the Barnett effect, and by Martin (1971), for magnetisation resulting from the 
grain's charge. 

Paramagnetic relaxation of the body in the external interstellar magnetic field had long 
been thought as the cause of the grain alignment. The "Barnett equivalent magnetic field", 
i.e. the equivalent external magnetic field that would cause the same magnetization of the 
grain material, is H BE = 5.6 x lCT 3 cj( 5 ) G, which is much larger than the interstellar magnetic 
field. Therefore the Barnett relaxation takes place on time scale tsar ~ 4 x 10 7 u;^ 2 , i.e., 
essentially instantly compared to the time of paramagnetic alignment ~ 10 11 s that establishes 
in respect to the external magnetic field 33 Thus Purcell (1979) concluded that J should be 
parallel to the axis of the grain major inertia. 

From the discussion above it looks that the precise time of internal dissipation is irrele- 
vant, if the internal alignment takes place faster than the alignment in respect to magnetic 
field. This is not true, however. Spitzer & McGlynn (1979), henceforth SM79 observed that 
adsorption of heavy elements on a grain should result in the resurfacing phenomenon that, 
e.g. should remove early sites of H 2 formation and create new ones. As the result, H 2 torques 
will occasionally change their direction and spin the grain down. SM79 showed that in the 
absence of random torques the spinning down grain will flip over preserving the direction of 
its original angular momentum. However, in the presence of random torques this direction 
will be altered with the maximal deviation inflicted over a short period of time just before 
and after the flip, i.e. during the time when the value of grain angular momentum is minimal. 
The actual value of angular momentum during this critical period depends on the ability of 
J to deviate from the axis of maximal inertia. SM79 observed that as the internal relaxation 
couples J with the axis of maximal inertia, it makes randomisation of grains during crossover 
nearly complete. Due to finite times of the internal relaxation the coupling is not complete 
and the grains preserve partially their alignment during crossovers. Calculations in SP79 
showed that the Barnett relaxation is fast enough to make randomisation of grains nearly 
complete. 

Lazarian & Draine (1997), henthforth LD97, revisited the problem of crossovers and 



32 A more exact statement is that the period of precession (about B ) of the magnetic moment M (and 
of Z aligned therewith) is much shorter than the mean time between two sequent flip-overs of a spinning 
granule (Purcell 1979, Roberge et al 1993) 

33 Radiative torques (Draine & Weingartner 1996, 1997) are currently favoured by observations as the 
mechanism of grain alignment (see Lazarian 2002). The alignment time for the radiative-torque mechanism 
is much larger than the internal dissipation time. 
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found out that the result is different if the thermal fluctuations associated with the relaxation 
mechanism are accounted for (see Lazarian & Roberge 1997). LD97 observed that the thermal 
fluctuations partially decouple J and the axis of maximal inertia and therefore the value of 
angular momentum at the moment of a flip is substantially larger than SM79 assumed. Thus 
the randomisation during a crossover is reduced and LD97 obtained a nearly perfect alignment 
for interstellar grains rotating suprathermally, provided that the grains were larger than a 
certain critical size a c . The latter size was found by equating the time of the crossover and 
the time of the internal dissipation t dis . For a < a c Lazarian & Draine (1999a) found new 
physical effects, which they termed "thermal flipping" and "thermal trapping" . The thermal 
flipping takes place as the time of the crossover becomes larger than t^ s . In this situation 
thermal fluctuations will enable flipovers. However, being random, thermal fluctuations are 
likely to produce not a single flipover, but multiple ones. As the grain flips back and forth, 
the regular (e.g. H 2 ) torques average out and the grain can spend a lot of time rotating 
with thermal velocity, i.e. being "thermally trapped". The paramagnetic alignment of grains 
rotating with thermal velocities is small (see above), and therefore grains with a < a c are 
expected to be aligned only marginally. The picture of preferential alignment of large grains, 
as we know, corresponds to the Serkowski law and therefore the real issue is to find the value 
of a c . In fact, Lazarian & Draine (1997) showed that a c ~ 1CT 5 cm for the Barnett relaxation 
mechanism. 

These finding stumulated more interest to the competing mechanism of internal energy 
dissipation, namely, to the inelastic relaxation of energy. Lazarian & Efroimsky (1999) redid 
the analysis of Purcell (1979) and obtained the expression for the relaxation time 

/ r \ 3/2 Oil q— 2.5 

* « «» |1 + (c/a)T (£) lf> kT^ , Q ^ (7-D 

where the grain is modelled by a 2a x 2c x 2c prizm (a < c ). Unlike our earlier 
expression (4-22), the kinetic energy of a rotating body is expressed using the factor of 
suprathermality, i.e. as (3kT gas . For crossovers thermal rotational velocities are important 
and for those Lazarian & Efroimsky concluded that inelastic relaxation may be dominant for 
grains larger than 1CT 5 cm, provided that the grains are substantially flat (e.g. axes ratio is 
4:1). The inelastic relaxation should also dominate atomic-size grains for which paramagnetic 
response fails (see discussion of paramagnetic relaxation within such grains in Lazarian & 
Draine 2000). As far as the critical-size value a c was concerned, the study by Lazarian & 
Efroimsky (1999) did not change its appreciably apart from the case of extremely oblate 
grains. 

However, in a recent paper, Lazarian & Draine (1999b) reported a new solid state effect 
that they termed "nuclear relaxation". This is an analog of Barnett relaxation effect that 
deals with nuclei. Similarly to unpaired electrons nuclei tend to get oriented in a rotating 
body. However the nuclear analog of "Barnett equivalent" magnetic field is much larger 
and Lazarian & Draine (1999) concluded that the nuclear relaxation can be a million times 
faster than the Barnett relaxation. If this is true a c becomes of the order 10~ 4 cm, which 
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means that the majority of interstellar grains undergo constant flipping and rotate essentially 
thermally in spite of the presence of uncompensated Purcell torques. For particles rotating 
at a thermal Brownian rate, nuclear relaxation dominates the inelastic one for chunks up to 
2 m. However, such large pieces rotate at a nonthermal rate and therefore inelastic relaxation 
rate that scales as u 3 compared with u 2 for the Barnett and nuclear relaxations is dominant. 

The scaling of inelastic relaxation with angular velocity ensures that this mechanism 
is dominant for fast rotating dust particles. Such a rotation is possible due to differential 
scattering of left and right hand polarized photons (Dolginov & Mytrophanov 1976, Draine 
& Weingartner 1996, 1997). This difference in scattering cross sections arises naturally for 
irregular grains (Draine & Weingartner 1996). Unlike the H 2 torques proposed by Purcell 
(1979), the radiative torques are not fixed in the body coordinates. Therefore they do not 
average out during fast thermal flipping of the grain. Hence, it is likely that they can provide 
means for suprathermal rotation for grains that are larger than the wavelength of the incoming 
radiation. The rate of internal relaxation for suprathermally rotating grains is important for 
special circumstances when the alignment happens very fast, e.g. in the stellar winds and 
comet atmospheres. 

Another interesting consequence of the coupling of rotational and vibrational degrees of 
freedom, resulting from the inelastic dissipation, is the expected relation between microwave 
and infrared polarisation arising from ultrasmall (a < 10~ 7 ) grains. Dipole emission from 
those grains was invoked by Draine & Lazarian (1998a,b) to explain the anomalous microwave 
emission of galactic origin within the range of 10-90 GHz. This emission was shown to 
interfere with the intended measurements of Cosmic Microwave Background (CMB), and 
therefore attracted a lot of attention (see Lazarian & Prunet 2002 for a review). In view of 
the attempts to measure the CMB polarisation, it is important to know whether or not the 
emission from small grains will be polarised. This problem has not had been fully resolved 
as yet. 

It is well known that the infrared emission from ultrasmall granules (the one of wave- 
length 12 fxm) takes place as the dust particles absorb UV photons. Can we check the align- 
ment of ultrasmall grains via infrared polarimetry? The answer to this question depends on 
the efficiency of internal relaxation. Indeed, UV photons raise the grain temperature and 
randomise grain axes in relation to its angular momentum vector (see Lazarian & Roberge 
1997) on the time scale of rotation- vibrational coupling. This is the time of internal relax- 
ation. Taking the necessary figures for Barnett relaxation from Lazarian & Draine (1999), 
we get the randomisation time of a 10~ 7 cm-sized grain to be 2 x 10~ 6 s, which is less than 
the grain cooling time. As a result, the emanating infrared radiation would be polarised 
very marginally. If, however, the Barnett relaxation is suppressed due to small-size effects 
(Lazarian & Draine 2000), the randomisation time is determined by inelastic relaxation and 
is ~ 0.1 s, which would entail a partial polarisation of the infrared emission 34 . The uncer- 



We would remind, that even though the inelastic relaxation of a dust granule is several orders more 
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tainty of this estimate arises from extrapolarion of dust properties and a classical treatment 
of an essentially quantum system. Further research in this direction is necessary. 

8. CONCLUSIONS 

1. In the article this far we have described the present situation in the studies of the 
dynamics of an unsupported top, and some of its applications to tumbling asteroids and 
comets and to the cosmic-dust particles. We addressed relaxation of excited (out of principal 
state) rotators through energy dissipation resulting from nut at ion- caused stresses. 

2. In many spin states of an unsupported rotator, dissipation at frequencies different 
from the nutation frequency makes a major input into the inelastic-relaxation process. These 
frequencies are overtones of some "basic" frequency, that is LOWER than the precession 
frequency. This is a very unusual example of nonlinearity: the principal frequency (precession 
rate) gives birth not only to higher frequencies but also to lower frequencies. 

3. In many spin states, the inelastic relaxation far more effective than believed hitherto. 

4- However, if the rotation states that are close to the separatrix on Fig.2, the lingering 
effect takes place: both precession and precession-damping processes slow down. Such states 
(especially those close to the homoclinic point) may mimic the principal rotation state. 

5. A finite resolution of radar-generated images puts a limit on our ability of recognising 
whether an object is nutating or not. Nutation-caused changes of the precession-cone half- 
angle may be observed. Our estimates show that we may be very close to observation of the 
relaxational dynamics of wobbling small Solar System bodies, dynamics that may say a lot 
about their structure and composition and also about their recent histories of impacts and 
tidal interactions. Monitoring of a wobbling comet during about a year after it leaves the 3 
AU zone will, most probably, enable us to register its precession relaxation. 

6. Measurements of the damping rate will provide us with valuable information on at- 
tenuation in small bodies, as well as on their recent histories of impacts and tidal interactions 

7. Since inelastic relaxation is far more effective than presumed earlier, the number of 
asteroids expected to wobble with a recognisable half-angle of the precession cone must be 
lower than expected. (We mean the predictions suggested in (Harris 1994).) Besides, some of 
the small bodies may be in the near-separatrix states: due to the afore mentioned lingering 
effect, these rotators may be "pretending" to be in a simple rotation state. 

8. Though the presently available theory predicts a much higher relaxation rate than 
believed previously, this high rate may still be not high enough to match the experimentally 



effective that was presumed in (Purcell 1979), its intensity (as compered to the Barnett relaxation) depends 
upon the granule size (Lazarian & Efroimsky 1999). For the smallest granules, the Barnett dissipation still 
dominates (Lazarian & Draine 2000). 
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available data. In the closemost vicinity of the principal spin state the relaxation rate must 
decrease and the rotator must demonstrate the "exponentially-slow finish". Asteroid 433 
Eros is a consolidated rotator whose Q-factor should not be too low. It is possible that this 
asteroid was disturbed sometimes in its recent history by the tidal forces. Nevertheless, it 
shows no visible residual precession. Hence, there may be a possibility that we shall have to 
seek even more effective mechanisms of relaxation. One such mechanism may be creep-caused 
deformation leading to a subsequent change of the position of the principal axes in the body. 

9. Inelastic coupling of the cosmic-dust grain's rotational and vibrational degrees of free- 
dom influences randomisation of grain axes when an ultrasmall grain absorbs a UV photon. 
As a result, the microwave dipole emission arising from ultrasmall grains may be polarised, 
while the near-infrared emission arising from the same grains may be unpolarised. 
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